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Abstract 

In this paper, we give a rigorous derivation of Einstein's geodesic hypothesis in general relativity. 
We use scaling stable solitons for nonlinear wave equations to approximate the test particle. Given 
a vacuum spacetime ([0, T] x R 3 ,/i), we consider the scalar field coupled Einstein equations. For all 
sufficiently small e and S < e q , q > 1, where 8, e are the amplitude and size of the particle, we show 
| the existence of solution ([0,T] x R 3 ,^,^) to the coupled Einstein equations with the property that 

the energy of the particle <f> e is concentrated along a timelike geodesic. Moreover, the gravitational 
field produced by e is negligibly small in C 1 , that is, the spacetime metric g is C 1 close to h. These 
results generalize those obtained by D. Stuart in [26], [27]. 

1 Introduction 

In general relativity, Einstein's geodesic hypothesis, which corresponds to Newton's first law of motion 
in classical mechanics, states that a free massive test particle will follow a timelike geodesic in the 
spacetime, where by free we mean in the absence of all external forces except gravitation, which is 
ascribed to the spacetime curvature instead of a force. For the concept of test particle, one has to ignore 
its internal structure as well as the gravitational field produced by it. This paper is devoted to a rigorous 
mathematical derivation of geodesic hypothesis in the sense that we approximate the test particle by 
real particles which are scaled stable solitons for a class of nonlinear Klein-Gordon equations. Since the 
real particles will interact with the background spacetime, we consider the following scalar field coupled 
Einstein equations 

\Rpv - \Rg^v = T^ig,^; V £ ,a), 

Here D g is the covariant wave operator for the unknown spacetime metric <? M „. 4> e is the complex scalar 
field representing the real particles. R^, R are the Ricci, scalar curvatures of the metric g respectively. 
V' e s is the first variation of the potential V e j. T^ v is the energy momentum tensor for the scalar field 4> e 
and is given as follows 

r^( ff) ^;Ve,5) =< d^,du4> e > -\g^ (< > + 2 Ve,5(<f )) , (2) 

where <, > denotes the inner product of two complex numbers, namely, < a,b >= |(a& + ab). 

The particle should have small size and small energy (or amplitude). We use the two small positive 
parameters e, 6 to denote its size and amplitude. In terms of local coordinate (t,x), the normalized scalar 
field (j> takes the following form 

<t>{t,x) = S~ 1 <p e (et,ex). 

Then the corresponding potential V for the normalized scalar field <fi is given by V e ^$ via the following 
equation 

V(0) = 6- 2 e 2 V e , 5 m- (3) 

In particular, one has V' e $((f> e ) — 5e~ 2 V '(S^ 1 ^). For a solution (Ai,g, <ff) of the above coupled Einstein 
equations, we find that the normalized scalar field cj> must satisfy the following nonlinear wave equation 

O gc( f, + V'(4>) = 0, g e (t,x)=g(et,ex). (4) 
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For small e, locally the above equation can be viewed as an equation on small perturbations of flat spaces. 
Since the particle will not change its shape during the motion, we use solitons of the above equation to 
approximate the normalized scalar field <\>. Note that the metric g € varies slowly, that is, \dg e \ < e||g||ci- 
By a soliton centered at point P in the spacetime, we mean a solitary solution of the form e lLOt f(x) to 
the equation n gC / P \<f> + V{4>) = on the flat space R 3+1 with metric g e (P). For concreteness, we will 
assume that the potential V takes the following form 

m = ^\<t>? - N(4>) - ^\<t>\ 2 - -^-M p+1 (5) 

2 2 p + 1 

with constants < to, f < p < 5. This type potentials guarantees the existence of solitons, which decays 
exponentially, see the work [1] of Berestycki-Lions. If <j> is close to such solitons, then the real particle 4> e 
is localized to a region of size e. Hence solitons centered at the position of the particles can be used to 
approximate the particles. 

We consider Einstein's geodesic hypothesis on a given vacuum spacetime which is diffeomorphic to 
([0,T] x R 3 ,/i) with the vacuum Einstein metric h. When the particle enters, it interacts with the given 
spacetime. To understand the motion of the particle, we consider the Cauchy problem for the scalar 
field coupled Einstein equations with initial data (R 3 , g, K, <j>\ ), where g is a Riemannian metric on 
R 3 , K is a symmetric two tensor, (<pQ,<fi\) are the initial data for the particle. The Cauchy problem is 
overdetermined which imposes compatibility conditions on the data: the constraint equations 

R(g) - \K\ 2 + {tvKf = |^| 2 + |V<^| 2 + 2V e ,«W§), V'X« - %trK =< V^g > . 

Here R(g) is the scalar curvature for the metric g and V is the covariant derivative with respect to g. 
We have the following main result 

Theorem 1. Let ([0,T] x R 3 ,/i) be a given vacuum spacetime for arbitrary T > with initial data 
(R 3 ,h, k). Assume 8 < e q , q > 1, 2 < p < |. Consider the Cauchy problem for the Einstein equations 
(1) with initial data (M. 3 ,g,K,cj)Q,<j)l) satisfying the constraint equations. Assume (g,K) is sufficiently 
close to (h,k) and the initial data (0q, 0f ) for the particle is close to some scaled stable soliton centered 
at the initial position Pq of the particle. Then for sufficiently small e, there exists a unique (up to 
diffeomorphism) solution ([0, T] x R 3 ,g, C ) such that <j> e is close some scaled soliton centered along a C 1 
curve which is close to a timelike geodesic starting from Po . Moreover, the spacetime metric g is close to 
the given vacuum metric h. 

Remark 1. Our proof will show that the same result still holds if S = for sufficiently small e$ which 
is independent of e. 

Remark 2. Existence of initial data satisfying conditions in the theorem will be shown in the last section. 

Remark 3. The lower bound of p (> 2) is required for regularity purpose: the spacetime metric has to 
be in C 1 . The upper bound p < | is needed to guarantee the existence of stable solitons. 

Remark 4. We have the same conclusion for much more general potentials V, for example, potentials 
V satisfying conditions given in [27]. 

A more precise statement of Theorem 1 as well as the definition of stable solitons can be found in 
the next two sections. We conclude from the theorem that in the limiting case when e goes to zero, 
the test particle moves along a timclike geodesic. This yields a rigorous derivation of Einstein's geodesic 
hypothesis in general relativity. 

Our theorem generalizes the result in [27] obtained by D. Stuart in two ways. First, our result 
describes the long time behavior of the test particle in the spacetime. In [27], it was shown that the 
solution (g,(j> e ) of the scalar field coupled Einstein equations exists only in a small portion [0,t*) x K 3 
of the given spacetime [0, T] x E 3 for some small positive constant t*. And the test particle travels 
along a timelike geodesic in the small time interval [0, t*). Here our result shows that we can extend 
the solution to the whole spacetime [0, T] x R 3 such that the energy of the particle <\f is concentrated 
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along a timelikc geodesic for arbitrarily large given time T. In particular, our theorem implies that the 
test particle moves along a timelike geodesic up to any given time T. Second, our result holds under the 
assumption S < e q , q > 1 which is weaker than <5 < e q , q > | imposed in [27]. Recall that 5 denotes the 
amplitude or the energy of the particle. If one wants to show that the gravitational field produced by 
the particle is negligibly small in C , i.e., \\g — h\\ c i < Ce, one has to show that the energy momentum 
tensor T^ v are bounded by Ce in £f s_1 , s > §. This requires S to be sufficiently small in terms of e. As 
the method for estimating the matter field (j) e € H s requires s to be an integer, D. Stuart proved the 
geodesic hypothesis under the strong condition q > |. On the other hand, notice that d^ 1 ^ is close to 
some soliton. In view of the equation (3), the potential V £i a, in particular the energy momentum tensor 
T^ u , has size 5 2 e~ 2 . Based on the heuristics that in the limiting case when e goes to zero, the potential 
V Cj 5 should be bounded, we see that the condition S < e q , q > 1 is needed to allow a proof of geodesic 
hypothesis in this setting. Our theorem thus answers the question of D. Stuart on the optimal value of q 
when q > 1. Furthermore, as was pointed out in Remark 1, the condition can be even weaker by merely 
assuming that S = eoe for some small constant e independent of e. This is robust and interesting as then 
the potential V t: s always has size e 2 , even when e goes to zero. 

The first aspect that we can extend the solution to arbitrarily large time T is based on a result on the 
orbital stability of solitons on small perturbations of Minkowski space. Note that the dynamics of the 
particle <\f are governed by the nonlinear wave equation (see the scalar field coupled Einstein equations 
(1)). In local coordinate (t,x), the corresponding normalized scalar field <p solves the nonlinear wave 
equation (4) on the slowly varying background with metric g e (t,x) = g(et,ex). The assumption on the 
initial data for the particle implies that the initial data for the normalized scalar field </> are close to some 
stable solitons. The property that the particle travels along a timelike geodesic can then be reduced to 
the orbital stability of stable solitons along a timelike geodesic on a slowly varying background. In other 
words, we need to show that the solution <fi of the nonlinear wave (4) is close to some soliton for all t and 
the center of the soliton propagates along a timelike geodesic. 

The related problem of stability of solitons in Minkowski space (g = m , the Minkowski metric), 
including the existence of solitons of the form e laJt / w (x) associated to the ground state of the corresponding 
elliptic equation 

Af OJ -(m 2 -u 2 )f u + \f OJ \P- 1 f u =0, |w|<m, 

has been studied extensively. For the ground state we summarize the relevant known results (see 

[1], [3], [19], [20], [24], [25] and references therein) in Theorem 2 in the next section. In particular, we have 
the C 4 regularity and exponential decay of the ground state fuj(x). The stable solitons are those with 
phase uj such that the associated energy ||| V/ W |j 2 2 ( R 3\ is convex in to, see [21]. Stable solitons have been 
proven to be orbitally stable for all time in Minkowski space first by J. Shatah [21] with radial symmetry, 
and later by M. Grillakis, J. Shatah and W. Strauss in [15], [16] in a much more general framework. More 
precisely, it was shown that if the initial data are close to some stable soliton in H 1 , then the solution 
to the nonlinear Klein-Gordon equation exists for all time and is close to some translated solitons(by the 
Lorentz transformation of Minkowski space). However, these works do not characterize the dynamics of 
the solitons. In particular, the centers of the solitons were not explicitly constructed. 

The first step along this direction was taken by M. Weinstein. In [28], he proved the orbital stability 
of solitons to nonlinear Schrodinger equations and gave the additional information on the position and 
speed of the solitons by using modulation theory. For the modulation theory, one decomposes the solution 
</> into the soliton part and the remainder part (j> — 4>s + v - The soliton part (f>s is unknown, depending 
on, e.g., its position and speed. Using the decomposition, the equations for 4> then lead to a linearized 
equation for the remainder v. We remark here that the linearized equation also depends on the unknown 
soliton (ps ■ We choose the decomposition such that the remainder part v is orthogonal to the generalized 
null space of the linearized equation at <f>s- This orthogonality condition together with the equations for 
4> leads to a coupled system of nonlinear ODE's (modulation equations), which governs the position and 
speed of the solitons. 
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Using this modulation approach, D. Stuart studied the stability of solitons to nonlinear wave equations. 
In [25] , he proved the orbital stability of stable solitons and showed that the center of the soliton moves 
along a C 1 curve which is close to a straight line in Minkowski space. Later in [26], D. Stuart studied 
the stability of stable solitons on small perturbations of Minkowski space. More precisely, he considered 
the Cauchy problem for the nonlinear wave (4) with initial data which are close to stable solitons on a 
slowly varying background with metric g e (t,x) — g(et,ex). He showed that stable solitons are orbitally 
stable and move along timelike geodesies up to time t*/e for some small positive constant t*. Due to 
the scaling, when applying this result to the problem of geodesic hypothesis, one can only show that 
the particle travels along a timelike geodesic in the short time interval [0,t*), sec [27]. Although t* is 
independent of e, it was required to be sufficiently small. 

A key improvement which allows us to obtain Theorem 1 is that we are able to extend D. Stuart's 
stability result up to time T/e for arbitrary large T > 0. We show that if the initial data are close to 
some stable solitons, then we can solve the nonlinear wave equation (4) up to time T/e and demonstrate 
that the solution is close to stable solitons centered along a timelike geodesic. The time T has to be fixed 
as we need to require e to be sufficiently small depending on T. However, we no longer need the smallness 
of T as it was in [26] . 

We will adapt the modulation approach to treat the orbital stability of stable solitons on a fixed 
slowly varying background. The new ingredient is that we can construct a coordinate system (Fermi 
coordinate system) such that the Christoffel symbols vanish along the trajectory of the center of the 
soliton, which is a timelike geodesic uniquely determined by the initial data. The motivation for choosing 
such a local coordinate system is to study the equation for the remainder part v by using the modulation 
approach mentioned above. Recall that we decompose the solution <j> of the equation (4) into soliton plus 
a remainder = <j>s + v. The remainder v is supposed to be small. To control v, we need to estimate 
n g e<j) S , in particular the term 



where is the Christoffel symbols for the metric g e . As the soliton <ps is expected to be centered along 
a timelike geodesic, to control the above term, a natural way is to choose a good local coordinate system 
such that along that geodesic the Christoffel symbols vanish. Furthermore, under such a coordinate 
system, the geodesic equations are linear. In particular, the geodesic can be parameterized by (t, u t) for 
some constant vector u e R 3 . This parametrization is exactly the one for straight lines in Minkowski 
space. Now as the full derivative of the metric components vanishes along the geodesic and the metric 
is slowly varying, that is, g £ (t, x) = g(et, ex), we conclude that near the geodesic, the metric g e is higher 
order (at least e 2 ) perturbation of the flat metric c/ £ (0,0). Hence the errors contributed by the soliton <j>s 
in the equation for the remainder v will have size e 2 . 

Nevertheless, for the orbital stability in the energy space H 1 on a slowly varying background, we are 
not going to study the equation for the remainder v directly. Instead, as in [25], we decompose the almost 
conserved energies of the full solution 4> around the soliton <ps- By using Gronwall's inequality, we can 
prove the orbital stability up to time T/e. The key is that we have avoided arguments based on boostrap. 
And by doing so we can remove the smallness assumption on t* which was used to close the boostrap 
assumption as it was in [26] . For the higher order Sobolev estimates of the remainder v which are needed 
to control the spacetime metric g, we turn to rely on the equation of v. 

We now briefly discuss the proof for the main theorem as well as the second aspect of our work that 
we can improve the amplitude S of the particle to be S < e q , q > 1 . We work under the Fermi coordinate 
system mentioned above. After scaling, it is equivalent to consider the scaled coupled Einstein equations 
for (g e , <fi) . Starting with the boostrap assumption on the unknown spacetime metric g e 
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V- det g' 



•s, 



\\d s (g*-h')\\ L 2 m (t)<2e 2 , l<s<3, t<T/e 



(6) 



4 



we can show that the normalized scalar field <j> is close to some solitons centered along the timclikc 
geodesic (t,u t) in H 3 . That is 



W(4>-<l>s)\\L>(R>)(t)<Ce, Vs<3, t < T/e 

for some constant C independent of e, where <f>s are solitons centered along (t, uot). Here we must clarify 
that we have used boostrap argument in order to estimate the spacetime metric g. But as discussed 
above, we avoid using boostrap argument when estimating the scalar field (p. The H 1 estimates are 
implied by the orbital stability of solitons. Then the higher Sobolev estimates follow by analyzing the 
nonlinear wave equation of the remainder <fi — 4>s- These Sobolev estimates for the scalar field <j> are used 
to control the the energy momentum tensor S 2 T flu (g e , </>; V(</>)) (after scaling) so that we can estimate the 
spacetime metric g e in order to close the above boostrap assumption. Choosing the relatively harmonic 
gauge condition [5], [17], one can turn the Einstein equations into a hyperbolic system for the components 
of the metric g e . Detailed reduction is carried out in Section 5. Combining with the vacuum Einstein 
equations for the given metric /i e , we can roughly obtain estimates for g e — h e as follows 

,T/e 

sup \\d a (f - h^\\ L2{m3) (t) < CS 2 / \\T^\\ H s(t)dt < CS 2 e- 1 

t<T/e Jo 

by using energy estimates for hyperbolic equations. This requires S < e q , q > | in order to close the 
boostrap assumption. Such a lower bound on q was suggested in [27] and the result there was proven 
under the even stronger condition q > |. 

However, through another new observation, we are able to improve the condition on S to be 5 < e q , 
q > 1 or 5 = e$e for sufficiently small eo which is independent of e. Notice that the tangent vector 
X = dt + uqOx of the geodesic (t, uot) is timclike and can be extended to a uniformly timelike vector field 
on the whole spacetime. We have already shown that the scalar field <j> decomposes into a soliton part 
4>s and an error term. Hence the energy momentum tensor 5 2 T flu (g e , <j>; V (</))) splits into the soliton part 
S 2 T^, which moves along the timelike geodesic (t,utot) or quantitatively 

8 2 \Xd s T^ u \<CS 2 e, Vs<3, 

and the error part 5 2 T^ V which satisfies the estimates 

5 2 \\d s T* || ia(R s)(t) < C5 2 e, s<3,t< T/e. 

To make use of the fact that the soliton part moves along the geodesic, when doing energy estimates, 
we multiply the hyperbolic equations by X(g e — h e ). Integrating by parts, we can pass the derivative X 
to the soliton part T^ v of the energy momentum tensor T^ u . As the soliton decays exponentially, using 
Hardy's inequality, we can show that 



/ 



T>„ ■ X(g< - h*)dx 

R 3 



< ||(1 + \x\)XT^\\ L2 \\(l + |x|)- V - V)\\ L . < Ce\W ~ & £ )IU*- 



Since the error part is small, the above observation allows us to improve the condition on 5. In fact, since 
the vector fields X, d t are uniformly timelike, the energy estimates for hyperbolic equations show that 

\\d(g e - h e )\\ 2 L2 (t) < CS 4 + CeS 2 f \\d{f - tf)\\ L .{ s )ds. 

Jo 

Here we have to require that initially \\d(g e — h e )\\ L 2(0) < CS 2 . Applying Gronwall's inequality or using 
boostrap argument, we obtain 

\\d(g*-h*)\\ L2 (t)<CS 2 , Vt<T/e. 
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Commuting the equations with d s , we have the same estimates for d s (g e — h e ), 1 < s < 3. Thus to close 
the boostrap assumption (6) and hence to conclude our main theorem, it suffices to require that 5 < e q , 
q > f or 5 = e e for sufficiently small e which is independent of e. Our main theorem then follows from 
the local existence result for Einstein equations in H 3 , see [4], [9]. 

Finally, we discuss the existence of the initial data (S , g, K,(f>Q, (/>{). As we have mentioned previously, 
we require the data to satisfy the constraint equations and the estimates 

\\d\g*-tf)\\ L ^ o) <C8 2 , Vl< s <3. 

For given <j>\), existence of initial data satisfying the constraint equations has been shown in [8], [5]. 
However, the above estimates do not follow directly from previous works. We will revisit the existence 
of initial data by using the implicit function theorem combined with the approach developed in [8], [6], 
[5]. To show that the data also satisfy the above estimates, we rely on a Hardy type inequality for a first 
order linear operator with trivial kernel in some weighted Sobolev space. We refer the reader to Lemma 
16 in the last section for details. 

In next section, we will address the basic setup, define solitons in Minkowski space and summarize 
the known results related to stability of solitons. In Section 3, we give the precise statement of our main 
theorem. In Section 4, we construct the good coordinate system such that the Christoffel symbols are 
vanishing along a given timelike geodesic. In Section 5, we reduce the Einstein equations to a hyperbolic 
system by choosing the relatively harmonic gauge condition. In Section 6, we prove the orbital stability 
of stable solitons along a timelike geodesic up to time T/e on a slowly varying background and show the 
higher Sobolev estimates for the solution which are used to control the energy momentum tensor T^ v . In 
Section 7, we prove our main theorem. In the last section, we discuss the existence of initial data. 

Acknowledgments The author is deeply indebted to his advisor Igor Rodnianski for suggesting this 
problem and for sharing numerous valuable thoughts, as well as tremendous indispensable help. 

2 Preliminaries and Stability Results in Minkowski Space 

In Minkowski space, for the nonlinear Klein-Gordon equation 

D^m^+r'^O, n = -d? + A = -d? + d 2 Xl +d 2 X2 +d 2 X3 , m>0, p>l (7) 

of complex functions <f>(t, x), looking for the solitons or stationary waves, that is, <j> is of the form e iU}t f ul (x), 
uj e R , we are led to consider the elliptic equation 

Af u -(m 2 -co 2 )f u + \f u f- 1 f u =0 (8) 

on R 3 . Such elliptic equation has been studied extensively and has infinite many solutions [2], of which 
we are particularly interested in the ground state, that is, solution of (8) with lowest energy 

Eu{v) = I f \Vv\ 2 + (m 2 - uj 2 )\v\ 2 - -l-\v\r +1 dx. 
2 J R 3 p+1 

It can be shown that the ground state has to be positive, radial symmetric. The existence of ground 
state as well as its properties are summarized in the following theorem, see [1], [3], [19], [20] , [24], [25] 
and reference therein. 

Theorem 2. For 1 < p < 5, uj G (— m, m), there exists a unique, positive, radial symmetric solution 
f u (x) e H 4 {M. 3 ) n C 4 (M 3 ) of the equation (8). It is decreasing in \x\ with the following properties: 

1 . Exponential decay up to fourth order derivatives 

|V%(a;)| < CV)e- c(aj)|a;| , VxeK 3 , \a\ < 4 (9) 
for some positive constant c(uj); 
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2 . Asymptotical behavior 

lim 1^- = -^^-, (10) 

3 . Scaling of the solutions 

f u (x) = (to 2 - L0 2 )^f(y/m 2 - lu 2 x), (11) 
where f(x) is the solution for m 2 — lo 2 = 1; 

4 ■ Identities of the energy 

3(p - l)(m 2 - lo 2 ) _ 5-p _ (p+ l)(m 2 - to 2 ) 

2iiv/.n 2 2 ~wm;~ u^ci, ■ ( } 

Existence of ground state has been obtained in [1], [3], [24]. K. McLcod [19] proved the uniqueness 
of the ground state. In [1], it was shown that f u (x) is C 2 and decays exponentially up to second order 
derivatives, which was generalized to be C 4 in [25]. The asymptotic behavior of the solution (10) was 
first proven in [20]. Using integration by parts, identities (12) follow by multiplying the equation (8) with 
U,x- Vf u respectively. 

Having the basic solitons e l0Jt f U! (x) of (7) corresponding to the ground state f u (x), one can study 
their stability. To start with, we must understand the symmetries of the equation (7), connecting to 
the symmetries of Minkowski space together with the structure of the equation (8), namely the scaling 
property of the ground state (11). These symmetries give an 8-parameter family of solitons. More 
precisely, denote A = (w, 9, £, u) g A with 

A = {(lj, u) g R 8 : |u| < 1, \ui\<m}. 

A subset of A is of particular importance 

A s tab = {(w,9,£,u) g A, |_^<^<i}, l<p< 7 -, (13) 

corresponding to the stable solitons. Define 

z(x; A) =A u (x-0=pP u (x-0 + (I-P u )(x-Z), p=(l-\u\ 2 )-^ (14) 
<d(x; A) = 9 — uu ■ z(x; A), 

where P u : M 3 — > M 3 is the projection operator in the direction u g R 3 , / is the identity map. Let 

<t> s {x;\)=j @(x ' X) U{z{x;\)), 

il> s (x; A) = e 10 ^ (ipuf u (z(x; A)) - pu ■ V z .Uz{x- A))) . 

Direct calculations show that <ps{x\\) solves (7) if the curve \{t) = (uj(t),9(t),t;(t),u(t)) obeys the 
evolution equations 

LO 

d) = 0, 9 = —, £ = u, ii = 0,. 
P 

Here we use the dot to denote the derivative with respect to t. We remark here that the centers (t, ut) of 
the solitons forms a straight line(geodesic) in Minkowski space. 
For A g A, let 

V(X) = (0,-,u,0). (15) 
P 

We find that V , s(A;x) = D\(f>s(\;x) ■ U(A)(inner product of vectors in R 8 ) and the important identity 

A x( f> s - m 2 s + |<fer Vs - D x ^ s ■ V{\) = 0. (16) 
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Here the Laplacian operator is for the variable x, used to distinguish the Laplacian A z for the variable 
z defined in (14). For a C 1 curve X(t) = (u(t),6(t),£(t),u(t)), let j(t) = (u(t),n(t),ri(t),u(t)) such that 

A = i + V(A). (17) 

We note that 



£(t) - / u(s)ds, n(t) = 9(t) - f 44 ds - 
Jo Jo P( s ) 



Consider the Cauchy problem for the equation (7) in Minkowski space with initial data (f>o{x) € 
_ff 1 (R 3 ), (f>i(x) e L 2 (R 3 ). The stability result for solitons is known. The following theorem is proven in 
[25]. 

Theorem 3. Let 1 < p < ^. For Ao G A sta b; there exists a small constant e(Ao) such that if 

e = ||0(O, a;) - Ao)|| ff i + \\d t (f>(0,x) - yj s (x; X )\\ L 2 < e(A ), 

then there exist a C 1 curve X(t) e A stab and a solution <f>{t, x) € H 1 of equation (7) satisfying 
\\<f>(t,x) -<p s (x;X(t))\\ H i + \\d t (f>(t,x)-i>s(x;X(t))\\ L 2 < Ce 

and 

\d t X(t) - V(X(t))\ < Ce 

for some constant C independent of e. 

Stability of solitons with Ao € A sta b has first been shown by J. Shatah in [21] for radial symmetric 
initial data and was later put into a very general framework in [15], [16]. Their approach relies on the 
fact that the energy E (jj (f ul ) is strictly convex in uj if initially Ao G A sta b, sec [21]. This condition on 
Ao is sharp in the sense that the solitons are unstable if the energy E u) (f UJ ) is concave in to, see [22], 
[23]. Alternatively, the modulation approach, pioneered by M. Weinstein [28] for showing the stability of 
solitons to nonlinear Schrodinger equations, leads to Theorem 3 which additionally gives the behavior of 
the curve X(t), see the work of D. Stuart [25]. 

We now briefly describe the modulation approach. Notice that equation (7) locally has a unique 
solution (cf)(t,x),dt4>(t,x)) € C 1 ([0, t*); H 1 x L 2 ). Decompose the solution <fi as follows 

<f>(t,x) = MvW)) + ^ x(t)) v(t,x), 

dt<t>(t, x) = fate A(t)) + e ie <* iA <*»u;(t, x) 
for a C 1 curve X(t) £ A sta b such that the following orthogonality condition hold 

< e^D^s, w > dx =< e~ ie D x ^s, v > dx , Vi e K. (18) 
Here in this paper for complex valued functions a(x), b(x), < a(x),b(x) >d x is short for 

ab + ab dx 



on IR 3 with measure dx. Differentiating (18) in t and using the equation (7), we can obtain a coupled 
system of ODE's for X(t). To estimate the curve X(t), we must control the radiation term (v,w). Wc 
first define two operators, L + and L_, appeared in the linearization in the real and imaginary part of 
the solution to (7). These two operators act on functions of z in i/ 1 (M 3 ), defined as follows 

i + = -A z + (m 2 - W 2 )- P /r 1 (z), 
L_ =-A z + (m 2 -uj 2 )-fP- 1 (z), 

which satisfy the following properties proven in [28] . 
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Proposition 1. We have 

(a) L- is a nonnegative self-adjoint operator in L 2 with null space kerL_ = span{/ w }; 

(b) L + is a self-adjoint operator in L 2 with null space kerL + = spa«{V z i/ w }| =1 . The strictly negative 
eigenspace of L + is one dimensional. 

It can be shown from the nonlinear wave equation of v(t, x) that the corresponding energy is 

E (t) = \\w + pu-\7 z v - ipujv\\ 2 L2(dz) + < Vi,L + vi > dz + < v 2 ,L-v 2 >dz, (20) 

where v, w are viewed as functions of (t, z). Although the operators L + , L- are not positive definite, one 
still can show that E (t) is equivalent to + IMU 2 under the orthogonality condition (18). 

Proposition 2. Assume A G A sta b- Assume v,w satisfy the orthogonality condition (18). Then there is 
a positive constant C(u,u), depending continuously on u, u such that 

C-V, U )(|H|| 2 + M\ 2 m ) < E Q (t) < C(u,u)(\\w\\ 2 L2 + \\v\\ 2 H1 ). 

We will use this proposition for our later argument. The proof could be found in [25], which is 
based on Proposition 1. The orthogonality condition < c~ ie Dgfis ,w >dx—< e~ ie Dgips, v >dx which is 
equivalent to < i(f>s,w > dx =< i4>Si v >dx shows that the energy E (t) is nonnegative. The proposition 
then follows by using a contradiction argument, see the detailed proof in [25]. Once we have control of 
+ ||HU 2 > by analyzing the ODEs for X(t), we can control the curve X(t) and obtain estimates for 
the solution </>(£, x). 

3 Statement of the Main Theorem 

Let (A4,h) be a vacuum spacetime. Assume M. is diffeomorphic to [0,T] x R 3 with coordinate system 
(t, x). We assume the vacuum Einstein metric h is C 4 and satisfies the following conditions 

^o _1/ 3x3 < (h kl )(t,x) < K I 3x3 , V(t,x) € [0,T] x R 3 , 

Kq 1 < -h 00 < K , \h^\<K , V/i,i/ = 0,l,2,3, (21) 
IWIc«([o,r|xR3) + \\\x\idh\\ L oo + IHarl^+^IUoo < K , \s\ < 2 

for some positive constant K , where h = h^ ly (t,x)dx fl dx u , x° = t, W{t,x) = (/i _1 ) AI!y (t, x), d is short 
for (d t , d Xl ,d X2 ,d X3 ). The first two lines require that the metric h is Lorentzian and the vector field d t is 
timclike. The decay assumption on dh is used to estimate g — h by using Hardy's inequality as initially 
we do not have estimates for \\g — ft,|| L 2( So )(our construction of initial data implies that for general data 
11.9 - ft|U 2 (Eo) is not bounded). 

The initial data (T,o,g,K,(j)Q,(j)l) for the Cauchy problem of system (1) with g a Ricmannian metric 
and K a symmetric two tensor on So, have to satisfy the following constraint equations on the initial 
hypersurface So 

R(g) - \K\ 2 + (trK) 2 = |^| 2 + |V^ | 2 + 2V £ , 5 (0 O ), m 
ViK ij -V i trK=<(l)\Si<i>l>, 

where V is the covariant derivative with respect to the Riemannian metric j on So- The potential V e .s is 
defined as in (3), (5). We remark here that <j>\ corresponds to n<\f for n the unit normal to S embedded 
to the solution (M,g). 

Since So is diffeomorphic to R 3 , we work under the coordinate system x on So- We assume the 
particle enters the spacetime at point £o <= S and is approximated by a scaled and translated stable 
soliton centered at £o- We hence has to define solitons on (M, h). The idea is that at any point P e A4, 
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wc simply define the solitons at P to be those on the flat space with metric h{P). More precisely, for any 
foliations [0,T*] x S T of M. with coordinate system (t,x), where S r C R 3 , let 

A st ab(t) = {(w,0,£,u h ) G K 8 | |— < £L < 1, ft +«fcflx is timelike at e [0,T*] x E T }. 
We define 

Definition 1. For a cwrue A(i) = (w(t), Ufc(t)) e A sta b(i); a soliton centered at in the 

direction Uh(t) on the space ([0, T*] x S T ,/i) is defined as follows 

<f>%(x; A(t)) = e . £ ->W-H f K(-(W)Q) /u(()(f -i 4g T (i _ 

w/iere 

m = a _1 (ai + UhQ), p = (1 - M 2 )~^, a > 0, a e R 3 , Q is 3 x 3 matrix 

We recall here that A u is 3 x 3 matrix defined in (14) gotc? m is the Minkowski metric. In particular, we 
can denote 

1> e s (x; A(i)) = ie^p^ua^six; A(t)) - u h {t)V x <t> e s {x; A(t)) 
corresponding to d t (j>% in the direction Uh and 

n<j) e s (x; X(t)) = ie _1 wp _1 ^s - uQ' 1 ■ V x <j) e s (x; X(t)) 
associated to n(f) e s , where n is the unit normal to the hypersurface S embedded to (M.,h) at time t. 

The above definition is well defined. We first show that h(t, £(t)) has a decomposition as in the 
definition, that is, a, a, Q exist. Notice that h is Lorentzian metric. At point the 3x3 matrix 

(h(t,£(t))ki) is symmetric and positive definite. We hence can find a 3 x 3 matrix Q such that 

QQ T = ((h(t,C(t)))ki). 

Then a, a are uniquely determined as follows 

a = (h 01 , h 2,h 03 )(t, £(t))(Q T )- 1 , a = y/-h(t, £(t))oo + aa T . 
Secondly, we prove that \u\ < 1. In fact, notice that the vector <9 t + Uhd x is timelike at (t, We have 

h(t, mho + 2h(t, £(t)) fc«£ + u k h h(t, £,{t))kiu l h < 0, 

which implies that 

\u\ 2 = a~ 2 (a + u h Q)(a + u h Q) T < 1. 

Finally, we demonstrate that the definition is independent of the choice of the 3x3 matrix Q. Let u, a, 
a, Q be another decomposition. Let P = Q~ X Q. We conclude that PP T — I and 

a = aP, a = a, u = uP. 

Hence we can show that 

u ■ (x - f (t))Q = uP-{x- £(t))QP = u-(x- £(t))Q, 
AuQ T = A uP P T Q T = P T A U Q T . 
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Since f u is spherical symmetric by Theorem 2, PP T — L, we thus have shown that (j)%{x; X(t)), ipg(x; X(t)) 
are well defined. 

For any function / on K 3 , we define the scaled Sobolev norm 

ll/lk = ||/(e)||*. = E el«l-f ||\7«/|U 2 . (23) 

|q|=0 

Additional to the constraint equations (22), we assume g, K, O , <j>\ satisfy the estimates 

||<TVoO«0 - 4> e s (x; \o)\\hi + eWS-^Kx) - ncf>%(x; X )\\ H , < C e, 

— (24) 
||V(fl - h)\\H ?( x o) + \\K k\\ H?(So) < CoS'e- 1 

for some constant Co and some point Ao = (coo, 9 , £ , u/ l (0)) <G A sta b(0), where h is the Lorentzian metric h 
restricted to the initial hypersurface So and k is the second fundamental form of So embedded to (Ai, h). 
In particular, Ao gives a timelike geodesic (t, 7o(t)) on (A4,h) such that 7o(0) = £o 5 7o(0) = u h(0)- 
Existence of such initial data set (So, g, K, <I>q,4>\) will be proven in the last section, see Theorem 6. 

Consider the Cauchy problem for system (1) with initial data (So,<7, K, <^) , (p{) described above. We 
have the following main result. 

Theorem 4. Let (Ai,h) be a vacuum spacetime satisfying (21). Assume 2 < p < |. Assume 5 < e g , 
q > 1 or 5 = e e. Suppose the initial data (T, ,g,K,(j)f > ,<j)\) satisfy conditions (22), (24). Then there 
exists e* > such that for all e, e G [0, e*) the Cauchy problem for system (1) admits a unique(up to 
diffeomorphism) solution (A4,g,(f) e ) with the following properties: there exists a foliation [0,T*] x S T of 
M with coordinates (s,y), where S T C R 3 and T* depends only on T, h and Ao, as well as a C 1 curve 
X(s) = (u}(s),0(s),£(s),Uh(s)) € A sta b(s) such that 

(1) The spacetime (A4 , g, £ ) is close to the given vacuum spacetime (A4,h) 

\\d{g-h)\\ Hl(i:s) <Ce, Vae[0,T*]; (25) 

(2) (j) e (t,x) is approximated by some translated soliton centered at £(s) 

||<rV e (s, y) - cj>%{y- A( S ))||h|(e s ) + («. v) ~ 4>Uv, < Ce, Vs e [0, T*]; (26) 

(3) The center of the particle £(s) is close to the given timelike geodesic (s,7 (s)) 

Us)- l0 {s)\ + \Lo{s)-LO \ + \u h {s)- 1 l Q {s)\<Ce, Vse[0,T*], (27) 

where the constant C is independent of e and the timelike geodesic (s,7o(s)) is the one starting from £o 
with speed u^(0). 

4 Construction of Fermi Coordinate System 

The geodesic hypothesis states that a test particle moves along a timelike geodesic. This timelike geodesic 
(i,7o(£)), using the notations in the previous section, is determined by the initial position £o and speed 
Uh(0). As having discussed in the introduction, our approach relies on the Fermi coordinate system such 
that the Christoffcl symbols are vanishing along the timelike geodesic (t,-f (t))- in this section, we give 
an explicit construction of such local coordinate system. 

Lemma 1. Let ([0,T] x M 3 ,ft.) be a smooth Lorentzian spacetime and (t, 7o(£)) be a timelike geodesic. 
Assume that the metric h satisfies (21). Then for all 8\ € (0, \T), there exists a subspace M 

[0, T - 2<5i] x R 3 c M c [0, T] x R 3 
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with a new coordinate system (s,y) <G [0, C(A'(0), h)(T — Si)} x M 3 such that the given timelike geodesic 
is represented as (s,u s) for some constant vector ii € R 3 , |uo| < 1- Moreover, along this geodesic, we 
have 

rjf>, u Q s) = 0, u Q s) = (mo) F1 Vs e [0, C(A'(0), ft)(T - «*„)], 

wftere T" are the associated Christoffel symbols and mo is the Minkowski metric, that is (too)oo = — 1 ; 
(wo)fcfe = 1, (ttiq)^ = 0,V/i 7^ za Furthermore, under the new coordinate system (s,y), there exists a 
positive constant K(5i 7 Jq(0)) depending on S\, 7 (0), ifo swcft iftai 

h(X,X)<-K(d 1 , 7 ' o (0)), X = d s +u k d Vk 

on the space [0, C(jq(0), h)(T — S\)] x R 3 . In particular, the vector field X = d s + u^d Vk is timelike. 

Proof. Let (t, x) be the given coordinate system on ([0, T] x R 3 , h). Our first step is to choose a coordinate 
system such that the given geodesic is represented by (t, 0). Let (t, x) be a new coordinate system defined 
as follows 

(t,x) = {t,x- 70 (*))• 

With this translation, we can show that the geodesic is parameterized by (t, 0). However, along this 
geodesic we only have r^t, 0) = 0. For simplicity, let's still use (t, x) to denote (t, x). 

Next, we change the coordinates in a small neighborhood of the geodesic such that all the Christoffel 
symbols are vanishing along the geodesic. To achieve this, we reproduce the construction here, which is 
essentially the same as that in [13], [18]. Let \{x) be a cutoff function 

X(x) = 1, |x| < r , 
X(x) = 0, \x\ > 2r 

for some positive number r . Introduce new coordinates defined as follows 

a* = x» + a£(t)x k X (x) + ^ kl (t)x k x l X (x), 

where (x°, x 1 , x 2 , x 3 ) — (t, x 1 , x 2 , x 3 ). Here we recall that Greek letters k,l run from 1 to 3 while the 
Latin letters fi, v run from to 3. The geodesic in this new coordinate system is still parameterized by 
(t, 0). The transformation law for Christoffel symbols 

~ _ dxP_dx^ ( dxT_ _ d 2 ^ \ 
m ~d&d&\ 01 dx^ dx^dx^J 

together with fact that ^Q (t, 0) = 0, implies that T£ v vanishes along the geodesic (t, 0) is equivalent to 
the following ODEs for a£(t), b^(t) 

/K = r o fc (*.o) + a^rj fc (t,o), 
\6& = r£(*.o) + ^rL(*.o)- 

Now take a^(0) = b kl (0) = initially. We conclude that the coefficients a k , b k[ are uniquely determined 
up to any time T. Moreover, we have the estimates 

K\ + \b^\<C(h,T, l0 ), 

where the constant C(ft,T, 70) depends only on the metric h, the time T and the geodesic (t, 7o(i)). 
Therefore for any 5i > 0, there exist ro sufficiently small and the associated cutoff function \ such that 

x : [0, T — 5i] x M 3 -> [0, T] x R 3 , i" 1 : [0, T - 2<5i] x R 3 [0, T - Si] x R 3 

are inclusions. Hence let M = x([0,T — <5i] x R 3 ) with the induced metric ft. We have 

[0, T - 26i] x R 3 c M c [0, T] x R 3 
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with a coordinate system x such that the Christoffel symbols vanish along the timelike geodesic (t, 0). 

We still need to change the coordinate such that ft — too along the geodesic. We have shown that 
under the coordinate system x = (t,x'), TP(t,0) = 0, which imply that dh(t,0) = 0, h(t,0) = ft(0, 0). 
Since the geodesic is timelike, we have ftoo(0,0) < and (ft(0,0))ju is positive definite. Assume 

ft(0,0)=( 'J a>0. 

Let Q3x3 be such that QPQ T = hx3- Consider the following new coordinate system 

(s, y) = (t, x>) (J aP j")( Va + £ f" aT Q °- x ) = (Wa + aP-^T, (taP^ + x')Q^). 

We can show that under this new coordinate system, the given timelike geodesic is parameterized by 
(s,u s). Moreover, h(s,u s) — too for some constant vector u € R 3 given as follows 

aP~ 1 0~ 1 
^o= = =£l 3 , KI<1- 

Va + a-P i cr 

By the previous construction, we have (s, y) € [0, (T — 5i)C(7 o (0), ft)] xR 3 , where we denote C(7 (0), /i) = 
Va + ttP _1 a T depending only on 7 (0) and the metric ft. 

Finally, notice that under the coordinate system x = (t,x'), the vector field d t is timelike, that is, by 
(21) and the construction of x 

h(d tl d t )<-K(5 ll7 ' (0)) 

for some constant depending on ft, Si, %{0). We remark here that the positive constant K(Si,%(0)) 
also depends on the cutoff function \- However, when 8\ is fixed, the cutoff function is also fixed. 
Therefore the corresponding vector field X = d s + u^d yk on the space (M, ft) with coordinate system 
(s, y) e [0, (T - Si Wa + aP^ 1 a T ] x R 3 is also timelike and satisfies 

h(X,X)<-K(S 1 , 7 ' (0)). 

□ 



5 Reduced Einstein Equations 

In this section, we reduce the scaled Einstein equations to a hyperbolic system for the components of the 
metric g € by fixing a relatively harmonic gauge condition, see [5], [17]. 

The Einstein equations are independent of the choice of local coordinate system. To solve the Einstein 
equations, we work under the Fermi coordinate system constructed in the previous section. More precisely, 
starting with the given vacuum spacetime ([0, T] x R 3 , ft) with local coordinate system (t,x), Lemma 1 
shows that we can find a new coordinate system (s, y) e [0, T*] x R 3 on M, where [0, T — 2Si] x R 3 c 
M C [0,T] x R 3 , such that the timelike geodesic is represented as (s,u n s) for some constant vector field 
u € R 3 , |t*o| < 1 and 

h(s, u a s) = too, T™„(s,u s) =0. 

Moreover, the vector field 

X = d S + U k Q dy k 

is uniformly timelike. To simplify the notations, we work on the space ([0, T] x R 3 , ft) for arbitrary T > 
with Fermi coordinate system (s,y). 

With this local coordinate system (s,y), to understand the particle 0% which has small size e and 
small energy 5, we instead consider the scaled Einstein equations. More precisely, recall that the potential 
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V e ,s satisfies the scaling V e ^(0) = S 2 e 2 V(S 1 <p). We infer that if (g(s,y),4> e (s,y)) solves system (1) on 
the space [0,T] x R 3 , then 

9 e (s, y) = g{es, ey), <j>(s, y) = <T V(es, ey) 
satisfy the rescalcd Einstein equations 

JiV(.9 e ) - ±R(g e )(g% v = 6 2 T^(g*, <j>- V (<(>)), 
\n fl .0 + V(0) = O 

on the rescaled space [0, T/e] x R 3 . Conversely, a solution (g e , <f)) of (28) on the space [0, T/e] x R 3 gives a 
solution (g, </> e ) of (1) on the space [0,T] x R 3 . It hence suffices to consider the scaled Einstein equations 
(28) on the space [0, T/e] x R 3 with initial data determined by (g,K,4>Q,<j)\) as well as the coordinate 
change from (t, x) to the Fermi coordinate system constructed in the previous section. 

Remark 5. The scaling heavily relies on the Fermi local coordinate system. However, the system (28) 
itself is independent of choice of local coordinate system on the space [0, T/e] x R 3 . 

We must determine the initial data for the corresponding Cauchy problem for the rescaled Einstein 
equations (28). By the construction of the Fermi coordinate system in Lemma 1 and by Definition 1, we 
define 

MV) = S-^Key), Mv)=8- 1 eft{ey), g e (y) = g(ey), K*(y) = eK(ey) 
on the initial hypersurface R 3 with coordinate system {y}. Denote 

\% = (cj , e^flo, 0, M ), n(t>s(y; Aft) = ipcj<t> s (y; Aft) - u Q- 1 V y <t>s{v\ Aft), 

where Ao = (coo, 9o, Co, «/i(0)) is given in Theorem 4 and Q is the 3x3 matrix in Definition 1. By (22), 
we have 

{Rm - |^| 2 + (trK^f = 5 2 (\^\ 2 + |Wo| 2 + 2V(0 O )), ( ) 

\ V'J% - V l trK" = S 2 < 0i, V^ >, 

where the covariant derivative is with respect to the metric g e . Using this scaling, the condition (24) 
becomes 



\\Mv) - My; WWh* + \\Mv) - ™My; K)\\h- < c e, 

\\V(g*-h*)\\ H2 + \\K*-k*\\ H2 < C S' 



(30) 



where h c {y) = h(ey), k e (y) = k(ey), V is the covariant derivative for the metric h e . We hence obtain 
the initial data (So, g e , K e , 4>o, 4>i) f° r rescaled Einstein equations (28). We also remark here that the 
existence of the initial data (£o,<7, K, (j>Q, <f>\) in our main theorem 4 is equivalent to the existence of 
(T, ,g e 7 K e ,(j)o,<j)i) satisfying the conditions (29), (30) under the Fermi coordinate system. Existence of 
such initial data set under certain conditions is discussed in details in the last section. 

The Cauchy problem for system (28) is underdetermined as any diffeomorphism of the spacetime 
([0,T/e] x M. 3 ,g e ,(f>) satisfying (28) leads to another solution. Such freedom can be removed by choosing 
a gauge condition such that system (28) is equivalent to a hyperbolic system for the components of the 
metric g e . Under the Fermi coordinate system on the space ([0,T/e] x R 3 , h e ), we define 

£ V, V) = (<? e r (i^, - f *„), a = o, i, 2, 3, 

where , are Christoffel symbols for the unknown metric g e and the given vacuum metric h e 
respectively. The gauge condition that we choose is 

9 x {g t ,h t )=0, VA = 0,1,2,3, (31) 
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which will be called the relatively harmonic gauge condition, sec [5], [17]. Instead of considering the full 
Einstein equations (28), we consider the following reduced Einstein equations 

R^if) - \R{9')9%, ~ \ (<£aV„0 A + «£ A V„0 A - 9^xQ X ) = 5 2 T^{g f -, 0; V(<j>)), (32) 

where the covariant derivative V is for the metric g e on the space [0, T/e] x M 3 . Then a solution of the full 
Einstein equations (28) can be constructed as follows: we write down (32) under the Fermi coordinate 
system and check that it is hyperbolic for the components of the metric g e . We then construct initial data 
(<? £ (0), dtg e (0)) from the given data (g e ,K e ) such that the gauge condition (31) holds initially. Hence we 
can get a unique short time solution ([0,t*) x R 3 , g e , </>) of the reduced Einstein equations (32) coupled 
with the matter field equation for (^(second equation in (28)). We then argue that the gauge condition is 
propagated, i.e., the relation (31) holds on [0,i*) x M 3 . This implies that the solution ([0,i*) x M 3 ,g e 1 (j>) 
to the reduced Einstein equations is, in fact, a solution of the full Einstein equations (28). 

Our first step is to check the hypcrbolicity of the reduced Einstein equations (32). We notice that 
under the fixed Fermi coordinate system(can be viewed as a local coordinate system), the full Einstein 
equations can be written as 

-(9T 9 a ^ v + gl x d, ((gT^ X a p) + 9^ ((sT^) + Q, u (f) = S 2 (2T,„ trT ■ g%), 

where 

Q,Ag e ) = c;^ v (f)d a g%d pg % x 

and C(g e ) denotes polynomials of g € a p, (</ e ) Q/3 . Taking trace of the reduced Einstein equations (32) , we 
get 

-R(g e ) + = trT. 

Plug this into (32). We have the reduced Einstein equations 

~(9T d aP g^ + gl x d» {{gT^U) + 9^ {(sT^p) + Q^tf) Q X d x g% = 5 2 {2T, V trT ■ g^). 
Similarly, we recall the vacuum Einstein equations for h e under the Fermi local coordinate system 

-(hT d a ^ v + K x d, (Vrr^) + h^ x d v (Vrr A „) + Q^h*) = o. 

Subtract the above two equations. We can show that the reduced Einstein equations coupled with the 
matter field equation for <f> are equivalent to the following hyperbolic system for the difference ip e — g e — h e 

-{9T P d a p% u + SP^ + SZ^ + SQ^ = S 2 (2T^ - trT ■ <£„), 
□ s .0 + V'(</>) = 0, 

where 

SQ^ = Q(s e ) - Qm - (gT^dxg^ + {h^t^d^, 

SZ^ = (fatf)"* - h^ihT") d v t^ + (gl x (gT - Kx(hT p ) - m af) d aP h^, (34) 
"V = ({gT^xg^ + W) Q/3 + faW) aP - m af> dxh%, - K x d^hT - ^ x d„(hT ) ■ 

Hence we have shown that the reduced Einstein equations (32) are hyperbolic. 

Secondly, we demonstrate that the gauge condition (31) propagates. Suppose (g e ,4>) satisfies system 
(33) on [0, t*) x M? for some small positive time t*. Take divergence of both sides of the reduced Einstein 
equations (32). Using Bianchi's identity and the fact that the energy momentum tensor T^ v is divergence 
free(due to the matter field equation of </>), we obtain the evolution equations for Q x 

V 1 V^ + [V iU ,V,]a' t = 0, i/ = 0,1, 2, 3, 
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where the commutator [V M ,V„] = V^V^ — V^V^. Hence Q x satisfies the above wave equations on 
[0,t*) x R 3 . Thus g x vanishes on [0,t*) x M 3 if Q», d t G^ vanish initially. We can make vanish on 
the initial hypersurface So by choosing initial data for g^, dtg^ (recall that only g e , K e or g\-, d t glj are 
given). Once we have — on E , we can show that d t Q^ also vanishes on E clue to the constraint 
equations. In fact, since (g e ,(j)) solves (32), the constraint equations (29) for g e , K e together with the 
vacuum constraint equations for h e , k e imply that 

V £„ + V„0o - 90,^^ = 0, v = 0, 1, 2, 3, 

where the covariant derivative V is for the metric g e . Hence d t G^ — on E if = initially. Therefore, 
we have shown that as long as (31) holds on E , a solution (tp e , (f) of the reduced Einstein equations (33) 
on [0,i*) x K 3 leads to a solution ([0,i*) x R 3 ,-0 e + h e ,(p) to the full Einstein equations (28). Since the 
solution of the Einstein equations (28) exists locally and is unique up to diffeomorphism [4], [9], it suffices 
to consider the reduced Einstein equations (33) with initial data (g e (0, y), d t g e (0, y)) such that (31) holds 
on the initial hypersurface So- 

It remains to construct initial data (g e (0), dtg e (0)) from (g e , K e ) on So such that the gauge condition 
(31) holds initially. Let 

(g e ) m (0) = -N- 2 , g e 0i (0) = Pu i = 1,2,3, 

where N is the lapse function and /3 is the shift vector on So- Since the Riemannian metric g e is the 
Lorentzian metric g e restricted to S and K e is the second fundamental form, we can show that 

gtj (0) = gtj , dtQtj (0) - ZNK^ + V^- + V,ft, 

where V is the covariant derivative for g e on So. That is (<7 e (0), dtgfj(0)) is uniquely determined by 
(g e ,K e ), N, fa. We claim that <9 t <7o M (0) is given by the gauge condition (31). In fact, the gauge condition 
implies that 

(2 - (mo)?)^ =2iV 2 (2(™o)£(gT - (^'(motf)^ + iV 2 (2(m )oV)^ - (sTVo)D^*^ 

+ N 2 g\,{gT^m Xv {2d a h^ - dvKp), r=g e - h e . 

We choose 9t<7§u on S as above. Notice that different lapse functions N and shift vectors (ii will lead to 
the same solution of the full Einstein equation (28) up to a change of local coordinate system. To better 
estimate the difference tp e , we simply set 

( 5 e ) 00 (0, y) = (h*) 00 (0, y), 5o ,(0, y) = h* oi (0, y). 

The above discussion shows that we have constructed the initial data (<? e (0), dtg e (0)) such that the gauge 
condition (31) holds initially. 

Finally, we show that g e is sufficiently close to h e (0) and hence is Lorentzian initially. Notice that 

dthtj = 2#(fcg)y + Vift + VjPi, 

where the covariant derivative is for the metric h e . By (21), (30), we can show that 

l|V( 5 e - /O||h 2 (e ) + \\W - ^)II^(E ) < C(h,C )5 2 , (35) 

where the constant C(h,Co) depends on h, Cq. Thus if e is sufficiently small, g e is a Lorentzian metric 
initially. We summarize what we have obtained in this section. 

Lemma 2. Let ([0,T] x R 3 ,h) be a vacuum spacetime. The Cauchy problem for (1) with initial data 
(So, g, K, cpQ, (j>\) is equivalentfup to a diffeomorphism and scaling) to the hyperbolic system (33), (34) 
with initial data (<po(0, y), </>i(0, y), g e (0, y), dtg e (0,y)) satisfying (30) and (35). 
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6 Stability of Stable Solitons on a Fixed Background 



We have shown in the previous section that to solve the full Einstein equations (1), it suffices to consider 
the hyperbolic system (33) on the scaled space [0,T/e] x R 3 with initial data (</>o(0, a;), </>i(0, a;), g e (0, x), 
dtg e (0,x)) satisfying (30) and (35). The standard local existence results imply that there is a unique 
short time solution (g e ,(j)) on [0,t*] x R 3 . Using continuity argument, the solution can be extended to 
T/e as long as (g e ,<j)) satisfies condition (30), (35) for a constant C independent of e. The estimate for 
g e — h e follows from energy estimates for hyperbolic equations if the source term T^„ — \trT ■ g € lies in 
C([0, T/e]; H 3 ). Since it is expected that the unknown metric g c is close to the given vacuum metric h e , 
the main difficulty for proving the main theorem is to show that (f> exists and is close to some translated 
soliton up to time T/e, that is, orbital stability of stable solitons along a timelike geodesic on the slowly 
varying background ([0, T/e] x R 3 ,g e ). 

In this section, we will let ([0, T] x R 3 , h) be a Lorentzian spacetime with the Fermi coordinate system 
(t, x) constructed in Lemma 1. More precisely, assume that (t, uot), t e [0, T] is a timelike geodesic, where 
the constant vector uq € R 3 , |t*o| < 1. Along this geodesic, we have 

h(t, uot) = mo, r° v (i, uot) = 0, 

where T^ v is the Christoffel symbols for h. Moreover, the vector field 

x = d t + 4d k 

is uniformly timelike on ([0,T] x R 3 ,/i). We assume the metric h G C 4 ([0,T] x R 3 ) and satisfies (21) for 
some constant Kq. Let g be another Lorentzian metric on [0, T] x R 3 . Consider the Cauchy problem for 
the rescaled nonlinear wave equation 

fn fl .0-mV + l0l p -V = O, (36) 
\<j>(0,x) = <j>o{x), d t <j>(0,x)=<j> 1 (x) 

on the rescaled space [0,T/e] x R 3 , where the slowly varying metric g e (t,x) = g(et,ex). We show the 
orbital stability of stable solitons to (36) if g — h vanishes as |x| — > oo and satisfies 

||d s+ V(V)llL 2(R 3) <2e 2 , |s|<2, Vte[0,T/e], (37) 

where ip = g — h, ip c — g e — h e . This condition can be viewed as a boostrap assumption. 
Our main result in this section is 

Theorem 5. Let g,h be Lorentzian metrics satisfying (21), (37) on the space [0,T] x R 3 with Fermi 
coordinate system (t,x). Assume (t,Uot) is a timelike geodesic for the metric h for a vector uq £ R 3 ; 
|«o| < 1- Assume 2 < p < |. Then for all Ao = (uio, do, 0, Uo) G A st ab; there exists a positive number e* 
depending on h, \ such that for all positive e < e* , if the initial data (f>o, 4>i satisfy 

\\<po(x) - <ps(x; A )|| ff i( R 3) + || 0i (a;) - ^s(x; A )||i,2( K 3) < e, (38) 

then there exists a unique solution <ft(t,x) € C([0, T/e]; H 1 (M 3 )) of the equation (36) with the following 
property: there is a C 1 curve A(t) = (co(t), 6(t), £(t) +u t,u(t) + Uo) € A sta b such that 

U(t,x)-<f>s(x;X(t))\\ Him + \\d t( l>(t,x)-M^m)\\L2 { m3) <Ce, Vte[0,T/e]. (39) 

Moreover 

|A(0) - A | < Ce, \j\ = \\-V(\)\<Ce 2 , Vie [0, T/e], (40) 
where j, V(X) are defined in (15), (17) and the constant C is independent of e. 

Remark 6. When g = h, D. Stuart in [26] proved the orbital stability result up to time t* /e for sufficiently 
small t* . 
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To show that the metric g e is C 1 , we need to control the energy momentum tensor T M „ in H 2 . We 
hence have to estimate the higher Sobolev norm of (f>. If initially </>o € H 3 , <f>i e H 2 , then we have 

Proposition 3. Assume 

ei = \\<Po( x ) - ^six; A )||h3 + \\(pi(x) - ip s {x; Xo)\\h 2 < oo. 
Let X(t) be the curve obtained in Theorem 5. Then 

\\d s (<j>(t,x) - fc&HtMLHW) <Cmax{e, ei }, Vte [0,T/e] 

kl<3 

/or a constant C independent of e, e\. 

A direct corollary of the above proposition is the boundedncss of the H 3 norm of the solution <j>. 
Corollary 1. If the initially data (j) <E H 3 , (f)\ <E H 2 , then 

\\d s Htr)\\m^)<C, Vte [0,T/e], 

|s|<3 

where the constant C is independent of e. 

The above boundedness of <fi in H 3 was used in [27]. In this paper, we have to use the fact that 
if initially the data (4>o,(f)i) are close to some soliton in H 3 , then the solution cf) stays close to some 
translated solitons for all t < T/e. 

Corollary 2. If ei < Cq£ for some constant Co, then 

]T W(4>-M*;Kt))\\voi>) < cc e, vte [o,T/ e ]. 

| S |<3 

To avoid too many constants, we make a convention that A < B means A < CB for some universal 
constant C depending on h, T, A , m, p. 

6.1 Decomposition of the Solution 

Since solution (f> of (36) exists locally, to begin with, we decompose the solution as follows 

U(t,o:) = 0s^;A(t))+e^( A ( t ))^^- W (t,x), d e (t,x) = (- d et g)i (et, ex), 
\d t( j ) (t,x) = Mx;X(t))+e^^ K ^w(t,x), p e (t,x) = y/-(g')°°(t,x), 

where A(t) € A. We also denote 



p e (t,x) = y/-(h e )°°(t,x), p(t,x) = ^/-h 00 (t,x). 

We define three functions depending on d € , p e , q e as follows 

ao = q € d € , ai=p f d £ , a = pj 1 q e = aoa^ 1 , b=p~ x d e . 

We now define the function q e (t, x) explicitly. Notice that h(t, uot) = mo and the Christoffel symbols 
for the metric h are vanishing along (t, uot). In particular we have dh(t,Uot) = 0. Hence we can show 
that 



\h" v (t,u Q t + x)\ < |(m )H +C Q \x\ 2 , V|x|<<5 , te[0,T], 

(42) 

?fc|", VFI ^ o , it [u, j j 

k=i 



\j2h kl (t,uot + x)^i > (l-Co\x\ 2 )J2\£ k \ 2 7 V\x\<5 , te[0,T] 
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for some positive constants Co, So, where we recall fi, v G {0,1,2,3} and k,l G {1,2,3}. Denote po = 
(1 — |uo| 2 )~ 5 . In particular, we have p > 1. Without loss of generality, assume \K^ 3 p^ 1 < 1 — 3poC'oS 2 , 
Ko > 10 and CoSqPo < jq- Then choose a function q(x) € C 3 (R 3 ) as follows 

!1 - 3C 0/ 5o|a;| 2 , M < ±6 , 
ix - 3 Po V(M 2 - f ) + (|«y - 2 - C a po)(5l \x\ 2 ), |x| e (|So, <J ), 

We define q e (t,x) as 

g e (i,a;) = g(e(a; - w i)). (43) 

We prove an inequality. 

Lemma 3. Let y = (2/1, y 2 , 2/3) € R 3 . T/ien 

m 2 C 2 fc 2 + qT^h^VkVl ~ 2(g £ p e )- 1 |(/ i e ) ofe 2/fc ||y • u\ - 3p m(q e p e )- 1 \(h e ) ok y k \\b\ > m 2 b 2 + \y\ 2 
for all iel, (t,x) G [0,T/e] xl 3 j = (2/1,2/2,2/3) € K 3 and u <= M 3 , |it| < 1. 
Proof. Fix t. After scaling, at point (/, uo^ + #), ^ suffices to show that 

mV 2 6 2 + q- 2 h kl y kyi - 2(qp)- 1 \h 0k y k \\y ■ u\ - Zpomiqp)- 1 \h 0k y k \ \b\ > m 2 b 2 + |y| 2 ,Vt G [0,T]. (44) 
When |a:| > <5 , we have q(x) — ^Kq Po ~ . Hence the left hand side of (44) 

> m 2 q- 2 b 2 + q^K^yf - 2V3q- 1 K 2 \y\ 2 - zVSpomq- 1 K 2 \y\\b\ 

> 15K$p 2 m 2 b 2 + (15Klp 2 - 8V3K 5 oPo )\y\ 2 - 12V3K*(%m\y\\b\ + m 2 b 2 + \y\ 2 
>m 2 b 2 + \y\ 2 . 

When \x\ < 5 , we have \h ok (t, u t + x)\ < C \x\ 2 by (42). Notice that q(x) < 1 - 3C poM 2 for \x\ < 5 . 
By (21), we can show that the left hand side of (44) 

> m 2 q- 2 b 2 + q- 2 (l - C \x\ 2 )\y\ 2 - 2q- 1 (l + C \x\ 2 )C \x\ 2 \y\ 2 

- SV^pomq-^l + Co|x| 2 )C |x| 2 |2/|H 

> (q- 2 - l)m 2 b 2 + q- 2 (l - q 2 - C \x\ 2 - 3qC \x\ 2 )\y\ 2 - Gpomq-^xf^ + m 2 b 2 + \y\ 2 

> 5C oP o\x\ 2 m 2 b 2 + q- 2 2C p \x\ 2 \y\ 2 - Spomq^C^x^yb] + m 2 b 2 + \y\ 2 
>m 2 b 2 + \y\ 2 , 

where we recall that Copo\x\ 2 < Co<5 2 po < jg- and p = (1 — |u | 2 )~^ > 1- Hence the lemma holds. □ 

In application, we need a similar inequality for the metric g e . 

Corollary 3. Let g e (t,x) — g(et,ex), h e (t,x) — h(et,ex). Assume \\g — h\\c° < e. Then 

m 2 q- 2 b 2 + q; 2 {gTvkVl - 2(q ePe )- 1 \(f) ok y k \\y ■ u\ - 2m(q ePe )- 1 \(g*) 0k y k \\b\ > (1 - Ce)(m 2 b 2 + \y\ 2 ) 

for all b G R, (t,x) G [0,T/e] x R 3 , y = (2/1,2/2,2/3) G R 3 and u G R 3 , \u\ < 1, where the constant C 
depends only on h and m. 
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6.2 Orthogonality Condition and Modulation Equations 

The decomposition of the solution (41) relics on X(t), which we write as 

X(t) = (w(i), 6{t), £(t) + wot, u(t) + uo). 
Denote 7(f) = (w(t),w(t),r](t),u(t)) such that 

A = 7 + U(A), 

where V(A) = (0, ^, u + uo, 0), p = (1 — \u + u \ 2 )~^ . Recall the notation defined in line (18). We choose 
X(t) such that the following orthogonality condition 

< D x fa, e i@ w > dx =< D X iP s , e m v > dx (45) 

holds. Differentiate it with respect to t. We can conclude that the orthogonality condition (45) holds if 
it holds initially and the curve A(i) satisfies 

< D 2 fa ■ \,e ie w > dx + < D x fa,d t (e m w) > dx =< D 2 ^ s ■ ke i@ v > dx + < D^ s ,d t (e l& v) > dx . 

Using the decomposition (41) and the relation A = 7 + V(A), we can show that the above equation is 
equivalent to 

(< D x ips,aoD X (t>s >dx - < D X (j)s,aiDxips >d x + < D 2 x fa,c ie w > dx - < D 2 x ip s ,e i@ v > dx )j 
= < D x (V(X)D X iP s ),e ie v > dx - < D x fa, oi(0« - D X ^ S V(X)) > dx (46) 
+ < D x ip s , a ((j) - fa) >d x - < aiD X (j) S + (ai - a )D x ip s , <j> t - ip s > dx , 

where we have replaced < D X V(X) ■ Dxfa,e te w > with < D X V(X) ■ D x ips,e l6 v > by the orthogonality 
condition (45). Denote 

H(t, x) = ba ge( f) + aid u 4> - A<j> = p^d^Ug^tj) + p e d e d tt (f> - A0 

= (pr 1 - (m f)d^ + 2 P ; 1 d e (gn ok d tk 4> + -^-d^dlifDd^ (47) 

where a^ k , are the corresponding coefficients. Recall the identity (16) for fa. Using integration by 
parts and observing that Re(e~ ie D x fa) = D x \fa\, we can show that 

< D x 4> s , A x 4> > dx =< D X (j> s , A x fa > dx + < A x D x fa, 4>-fa >d x 

=< D x fa, m 2 fa - Ifaf^fa + D x fa ■ V(X) > dx + < D x A x fa,<j> - fa > dx 
=< D x (D x ip s ■ V{X)),<j) - fa >dx + < D x fa, D x fa ■ V(X) > dx 
- < Dxfa, -m 2 + \<j>s\ p - 1 <t> + (p - l)\<Ps\ p - 2 <j>sRe{c- tB {4> - fa)) > dx ■ 

Now use the identity (47) to replace p e d € d t t<t> in (46). The equation (36) of <j) then implies that the right 
hand side of (46) can be written as 

F(t; A(i)) = < d D x ip s , <t> - 4>S >d x + < (a - a{]Dx^s - ^Dxfa, <j> t - fa >d x 

+ < D X (V(X)D X ^ S ), (oo - l){4> - fa) >dx + < D x fa, { ai - l)D x fa ■ V(X) > dx 

+ < d k (a k »D x <p s ) - WDxfadpt > dx - < bDxfa,c l& Af(X) > dx (48) 

+ < (b - l)Dxfa,m 2 ^ - |0s| p_ V > dx -<(b- l)faDx\fa\ p -\<f> - 4>s > dx , 

where we denote 

Af(\) = e- je (|^rV - Ifa^fa - \4>s\ p -\4> - <f>s) - (p - l)|^| p - 1 e ie i?e(c- je (0 - fa))) (49) 
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as the nonlinearity depending also on A. Here recall that b — p e 1 d e . In particular if v = v\ + iv 2 for real 
functions v\, v 2 , then Re(e~ ie ((j) — <f>s)) = (qed e )~ 1 vi by the decomposition (41). 

To further simplify the above modulation equations, denote 

D =< D x ^s,D x <t>s >dx - < D\<ps,D\ip s >dx, 

Dx =< L>aVs, (ao - l)D x <t>s >dx - < D x <j>s, (ai - l)D x ^s >dx, (50) 
D 2 =< Dl(f> s ,e ie w > dx - < DliP s ,e m v > dx . 

We point out here that D, Di, D 2 are 8x8 matrices. Hence we choose the curve X(t) such that 

(D + D 1 +D 2 )j = F(t;X(t)), X = j + V(X), 7 (0) = A(0) (51) 

and we require A(0) satisfies the orthogonality condition (45). 

Equations (51) are called the modulation equations for the curve X(t). If the orthogonality condition 
(45) holds initially and X(t) solves the ODE (51), then (45) holds as long as <f), X(t) exist. 

6.3 Initial Data 

We have reduced the orthogonality condition (45) to a coupled system of ODE's for A(t) if initially A(0) 
satisfies (45). We use the implicit functional theorem to show the existence of the initial data A(0) 
satisfying the orthogonality condition. 

Lemma 4. Denote 

e = \\M X ) - 4>s{x; A )||hi( R 3) + ll^i (x) - Ao)||l2 (r3) 

for some Ao € A s t a b- Then there exists a positive constant ei(Ao), depending only on Ao, such that if 
e < ei(Ao), then there exists A(0) € A sta b with the property that if 

Mx) = <t>s{x; A(0)) + e i9W0)) ^«(x; A(0)), 

M*) = Mx;X(o)) + e^^^^wix-Aio)), 

then the orthogonality condition holds 

< D^ s {x- A(0)), e ie ^v(x; A(0)) > dx =< D x <t> s {x; A(0)), e i& ^w(x; A(0)) > dx . 
Moreover, we have 

|A(0) - A | < C(A )e, 

\\v(x;X(0))\\ H i + ||t£;(ar; A(0))|| ia < C(X )e 

for some constant C(Ao) depending only on Ao- 

Proof. Define a functional T : H 1 x L 2 x M 8 ->• M 8 such that 

T(v, w, X) = < a 1 D x <ps(X;x),(j) 1 (x) - ips(X;x) > dx - < a a D x ips(^;x),(f) a (x) - (f>s(X;x) >dx ■ 

In particular, we have J"(0, 0, Ao) = 0. Notice that 

F x (0, 0, A ) =< a D x ip s , D X (f> s > dx - < ai D X (f> s , D x ^ s > dx = D + Di. 

By Lemma 6 and Lemma 5 proven later, we can conclude that if e is sufficiently small, then D + D\ 
is nondegenerate initially. Since F is Lipschitz continuous in (v,w), the implicit function theorem then 
implies that there exists A(0) G A sta b satisfying the orthogonality condition (45) as well as the estimates 
in the lemma. □ 
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6.4 Boostrap Argument 

By Lemma 4, we can choose A(0) close to Ao such that the orthogonality condition (45) holds initially. The 
local existence result shows that there is a short time solution <p(t,x) of (36). To prove the existence of 
solution A(i) of the modulation equations (5 1 ) , we have to demonstrate that the 8x8 matrix D + D x + D 2 is 
nondegenerate, which depends on X(t) itself. We thus need to consider the modulation equations coupled 
to the nonlinear wave equation (36). Since initially the radiation term (v,w) is small in H 1 x L 2 , we 
show that for the coupled equations, (v, w) stays small in H 1 x L 2 for all t <T/e, which implies that the 
modulation equations are solvable and we can control the modulation curve X(t). 

Proposition 4. Let (4>(t, x), X(t)) be solutions of (36), (51) on [0,T/e] x R 3 . Assume the complex 
functions v(t,x), w(t,x) satisfy the decomposition (41). Then for sufficiently small e, we have 

IK^aOlljn^s) + \\w(t,x)\\ L 2 {R3) < Ce, Vt G [0,T/e], 
\i(t)\<Ce 2 , Vt€[0,T/e] 

for some constant C independent of e. 

We use boostrap argument to prove this proposition. To estimate | | , we rely on the modulation 
equations (51) together with the fact that D + D\ + D2, as an 8 x 8 matrix, is nondegenerate. It turns 
out that D depends only on to and is nondegenerate for u> = luq. Thus if ui(t) is close to ujq, we have the 
nondegeneracy of D. From this point of view, we take a subset of A sta b defined as follows 

p — 1 ui 2 

A <5 = {{u,0,£,u)\w £ [w - S ,uj + 5 ] C {uj\ < — < 1}, \u - u \ < 5 }, (52) 

d Ap in 

where Sq is a positive constant, which will be determined in Lemma 9. Here we recall that A = 
(u>o, Oq, 0, u ) e A sta b. Thus we can choose So sufficiently small such that A^ is nonempty. Our first 
boostrap assumption is that X(t) e A,5 for all t £ [0, T/e]. 

By the definition of D\, D2 in line (50), we see that D\ has size e if the center (i, £(t) + u^t) of 
the soliton <f>s(x; A (t)) docs not diverge far away from the timelike geodesic (t,u t). Hence we expect 
that |£(t)| is uniformly bounded, which is also suggested by Proposition 4. In fact since £ = u(t) + r), if 
I7I = \(oj, 7T, r], u)\ < Ce 2 , then 

\m\ < \m\ + Ho)\t + ct 2 e 2 < \m\ + (c + i)t 2 , vt < r/ e . 

We remark here that this is compatible with Theorem 4 as if we scale it back to the space [0,T] x M 3 , 
the center of the soliton becomes (t, e£ + uot) which is close to the geodesic (t, u t). Similarly, D 2 is an 
error term if IHI^rs) + ||u|| L 2 (R 3 ) is small. 

To prove Proposition 4, in addition to the assumption that X(t) £ As , we assume 

< 2C 2 , Vte[0,T/e], (53) 
\\w(t,x)\\ L 2 {R3) + \\v(t,x)\\ H i {R3) < S u Vt £ [0,T/e] (54) 

for some constants C2, Si which will be fixed later on. Without loss of generality, we assume C2 > 1, 
S\ < 1, C|e < 1. These are the boostrap assumptions in this subsection. 

As having mentioned previously, we do not have estimates for \\g e — /i c ||l 2 (s ) initially. In fact, our 
construction of initial data implies that g e — h e is not bounded in L 2 for general data. To bound g e — h e , 
we rely on Hardy's inequality. 

Lemma 5. Let f(x) £ C 1 (M 3 ). Assume f{x) — >• as \x\ -4- 00. Then 

Wfix^W^^ < 6\\Wf(x)\\ mm . 

This inequality can be proven by using integration by parts under polar coordinate system. Detailed 
proof could be found in [12], [29]. In particular, using Sobolev embedding, the assumption (37) implies 
that 

II ^Vi(RS) + H^ £ (l + N)" 1 H^( R 3) + ||9 s+ Vll^( R 3)<e 2 , |s|<2, Vte[0,T/e]. (55) 
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To prove Proposition 4, we use boostrap argument. We show that the matrix D + D\ + D 2 are 
nondegenerate initially. Thus the modulation equations can be solved locally. If under the boostrap 
assumptions (53), (54), we can show that Proposition 4 holds for some constant C which is independent 
of C 2 , then we can close the boostrap assumptions if we take C2 = |£(0)| + (C + 1)T 2 , 5\ = 2e, <5 = 2CTe. 
We thus can conclude Proposition 4 if e is sufficiently small. 

The strategy is as follows: we first show the nondegeneracy the modulation equations and obtain the 
estimates for \ j(t)\. Under the orthogonality condition (45), we use energy estimate to demonstrate that 
the radiation (v,w) is small in H 1 x L 2 . 

6.4.1 Nondegeneracy of the Modulation Equations 

To obtain estimates of X(t) and to show the local existence of the modulation equations (51), we demon- 
strate that, under the boostrap assumption X(t) e A«5 , the leading coefficient D in (51) is nondegenerate. 

Lemma 6. Let D be the 8 x 8 matrix defined in (50). If X = (oj,9,£,u) € As , then 

\detD\ > C{S , A ) > (56) 

for some constant C(So, Ao) depending on 60, Xq. In particular, D is nondegenerate. 

Proof. Direct calculations show that 

D u9 = D u {cj\\f u \\ 2 L2 ), D u z = puD u B, D^ u = -pB(I + p 2 u-u), 
Duu = Dq£ = = D uu = D 6u = D eg = 0, 

where we denote 



B(t) - uj 2 \\Ux)\\ 2 L2{m + -||V./ w (z)||| 2(K3) . (57) 



1 

g II v xJu>\^}\\L 2 (BL3) 

Notice that D is antisymmetric. We can show that 



detD = |A,M|/ W ||!2)| 2 • |det^ 5u | 2 = |£) w (w||/ w ||| 2 )| 2 BV°- 
Using the scaling property (11) of f u (x), when < ^ < 1, we have 

D u (u>\\f u \\h) = (™ 2 ^y^ 2 )(™ 2 - ^ 2 )^ll/lli 2(R 3 ) < 0. 

Since p > 1, the lemma then follows. □ 
6.4.2 Estimates for the Modulation Curve 

We have shown that D is nondegenerate. To estimate ^(t) by using the modulation equations, we 
have to show that D\, D2, F(t;X(t)) are error terms. Estimates for D\, D2 can be obtained by using 
Cauchy- Schwartz inequality. The main difficulty for estimating F(t;X(t)) lies in the nonlinear term 
< D\(f)s,e ie Ar(X) >dx- Before touching the equations, we prove two lemmas. The first lemma gives 
control of the nonlinearity A/"(A). Let 

Nix) = — 
V 1 p+V 1 

for complex number x. In particular, we have D S N = i|x| p_1 a;, where x is the complex conjugate of x. 
By the definition (49) of the nonlinearity A/"(A), wc can write 

M{X) = 2D s N(f ul + (qed^v) - 2D s N(f LO ) - 2D i DN{f u ) ■ {q.d^v, 

where DN ■ v = D x Nv + D s Nv. 
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Lemma 7. Let Af(X) be defined in line (49). Assume <j> decomposes as (41). For all p > 2, we have 

WW\<\v\ 2 + \v\ P , (58) 

<M 3 + M P+1 - (59) 



N(f u +v)- N(f u ) - DN(f u )v - ^vD 2 N(f u )v 



Proof. Let v = (q e d £ ) x w. We can write Af(X) as an integral 

7V(A) = 4 f f sv 2 D 2 D s N(f ul + tsv)dtds. 
Jo Jo 

Since p > 2, we can show that 

\D 2 D s N(f u + tsv)\ < 1 + \f u + tsvf- 2 < 1 + \v\v- 2 . 
Hence (58) holds. The second inequality (59) follows similarly. □ 

This lemma implies that the nonlinearity in F(t; X(t)) is in fact higher order error term. For the other 
terms in F(t; A(t)), observe that if g = h, then p e — l,d € — l,q e — l have size e 2 near the geodesic (t, uot) 
by (42) and the fact that <f>s decays exponentially. To pass h e to g e satisfying (55), we use the following 
lemma. 

Lemma 8. Let F\{x), F 2 (x) be two C 1 functions such that Fi(0) = F 2 (0) = 0. Then we have 

WF^six; A))(W) - F 2 (m ))\\ Lr{R3) < H^iH^ ||F 2 || C1 (1 + \tf)e 2 , Vr e [1,2], 

where A = (cj, 6, £ + u t, u + u ) e A stab . Similarly 

WFtiMvWftV&WLrp*) < \\Fi\\d \\F2\\ C i(l + lCl)e 2 , Vr e [l,cx>]. 

Proof. Recall that <j>s{x; A) = c te f LU (z), z = A u+Uo (x — £ — u t). We conclude that |z| 2 > \x — £ — uot\ 2 . 
Hence by Theorem 2, we have 

|fi(<Mz;A))| < H^Hcie-^l'l < H^H^c-^-^l. 

Since dh(t, uot) = 0, h(t, uot) = mo, we can show that 

\F 2 (g')(t,u t + x) -F 2 (m )(t,u t + x)\ < \\F 2 \\ cl (e 2 \x\ 2 \\h\\ C 2 + \g e - h e \). 

Then (55) and Holder's inequality imply that 

WF^x; X))(F 2 (g<) - F 2 (m ))\\ Lr{R3) 



< + 1*1)11^ ||F 2 || C1 1|^(1 + \z\)-'\\ L , + ||^|| C1 ||F 2 || C2 (1 + ICIV 



lie 1 ii'^ i 1 1- iii 2 f ii- r iiic i ii- r 2ncH- L 1- i?' |2v ' 

< 11^11^11^11^(1 + |£| 2 )e 2 , re [1,2]. 

Similarly, by (55), we have 

WF^six; A))F 2 (^)|| ir(M3) < \\Fi(4>s{x; X))F 2 (dh e )\\ L r( R 3) + ||*2|MI 

<||^|| cl ||^|| cl (l + |e|)e 2 , re [l,oo]. 



□ 



Having proven the above two lemmas, we are able to estimate A(t). It suffices to show that £>i, Z) 2 , 
F are error terms. 
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Proposition 5. Under the boostrap assumptions (53), (54), we have 

Pill <(l + l£lV<(C 2 £) 2 , (60) 

IIAjII < HU»(r3) + NI^(r») < Ji, (6i) 

ll^ll<(C 2 e) 2 + || V ||^. (62) 

If e, <5i is sufficiently small, by Lemma 6, estimates (60), (61) show that D dominates D\ + D 2 . Hence 
D + Z?i + L> 2 is nondegenerate. Then using (62), we can estimate 7. 

Proof. Inequality (60) follows from Lemma 8 and the boostrap assumption (53). Estimate (61) for D 2 
can be obtained by using Cauchy-Schwartz's inequality and the assumption (54). 

For (62), we apply Lemma 7 to control the nonlinearity < bD\(j)s,e ie J\f(X) >d x and use Lemma 8 to 
estimate other terms by observing that 

b-l=p~ 1 d e -l, a - 1 = q e d c - 1, ai-l=p e d e -l, a k », b^ 

can be written as the form Fi(g e ) — Fi(mo). Hence we can show that 

\\F\\ < 6 2 (1 + ICI 2 ) + IIM 2 + \v\>\\ L i < {C 2 ef + \\vf Hl . 

□ 

To make D + D x + D 2 to be nondegenerate, we choose Si in the following way. Fix Sq. Then let e, Si 
be sufficiently small such that 

||I>i|| + ||I>2|| <C{C 2 ef+C5 1 < ^C(S ,X )^ 

where C is the implicit constant in Proposition 5, which by our notations is independent of e, C 2 . Here 
C(5 ,X ) is the constant in Lemma 6. So far, only <5o, C 2 arc unknown constants. However, all the 
implicit constants are independent of C 2 , e. With this choice of Si, we can estimate 7. 

Corollary 4. Let 5i be chosen as above. Suppose X(t) G A$ and^(t), w,v satisfy the boostrap assump- 
tions (53), (54). Then we have 

\i\<\\F\\<(C 2 ef + \\v\\ 2 Hl . 

This corollary shows that as long as the radiation term v in the decomposition (41) of the solution 
(41) is small, we can solve the modulation equations (51) and obtain estimates for the modulation curve 
j(t). The modulation equations are used to guarantee the orthogonality condition (45). Next, we show 
that under the orthogonality condition, the energy \\v\\ H i + \\w\\ L 2 of the radiation term (v,w) is small. 

6.4.3 Energy Decomposition 

We use energy estimates to show that the radiation term (v, w) is small in H 1 xL 2 . We consider the almost 
conserved energies for the full solution (f> of the nonlinear wave equation (36). Using the decomposition 
(41), we decompose the energies for <p around the solitons 0s- The associated energy for the solitons <f>s 
can be computed explicitly up to an error. Then combining with the orthogonality condition, we can 
obtain estimates for + |M|l 2 - We first define almost conserved energies for 4> and decompose them 

around the solitons <ps- 

We recall the energy momentum tensor with respect to the metric g e 

T M „[0] =< > ~g%{< 0^,0^ > +2V(0)), 

where 
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For a vector field Y and a real function f3, we have the identity 

D»(f3T llv [<f>}Y») = (3T>"'l<f>}n^ + f3< □^-mV+ICVW) >+T(Y,D/3), (63) 

where tt^ v = \dyg i tiV is the deformation tensor of Y and -D/3 is the gradient of the function (3. Consider 
the region R 3 x [0,t\. First we take (3 = 1, F = 8t- Using Stoke's formula and equation (36), we have 

H{t)=H{0)- f I T^[^TT%dvo\, (64) 
Jo Jr 3 

where 

H(t) = - f -T 0v [cj)]{g"f v da = \ [ < 8 k ^8 k <P >-< 8*0,8^ > +2V{<P)da. (65) 

JR 3 2 JR 3 

Here we recall that da = ^— det g € — d 2 dx, dvol — d 2 dtdx. 
Then let [3 = a = p^ 1 q e , Y — 8 k . We have 

n k (t) = n fe (0) + /* / aT^[<j>]n% + d^aTMdvol, (66) 
Jo Jr 3 

where 

n fc (t)=/ aT lik [d>](g e ) ^da= [ p- 1 q £ < 8*0,8^ > da. (67) 

JR 3 JR 3 

We remark here that life are conserved quantities of the equation (36) if the metric g e is flat. For general 
slowly varying metric g e , we will show in the next section that these quantities are almost conserved, 
that is, the error is of higher order. 

Since (j> is complex valued, we define the charge of the solution <f> 

Q(t) = / a < id l (t), 4> > da, da = d\dx. (68) 

JR 3 

Using the equation (36), we can get an integral form of Q(t) 

Q{t) = Q(0) + / / |(ac^) < <j) > +ad 2 (< idtd^, <p > + < i&fa 8 t <p >) dxds 
Jo Jr 3 8t 

= Q(0) + [ [ d\ < id^ad^fa 4> > +ad 2 < id^^, 8^ > dxds (69) 

JO JR 3 

= Q(0) + I I < id^ad^, <p > d 2 dxds, 

Jo JR 3 

where we have used the equation 

8^4> + d- 2 8„(d 2 e )d^4> - m 2 4> + \<j>\P- V = 
together with the fact that the quadratic form < i(g e )-, ■ > is antisymmetric. 

Remark 7. The reason that we put some weight in the decomposition (41), also in the definition of the 
almost conserved quantities Tl k , Q, is to reduce the positivity of the energy of the radiation term (v,w) 
to the case in Minkowski space, which has been proven in [25], see Proposition 2. 

Next, we expand H, H k , Q, as functionals of the full solution 0, around the soliton <f>g. The soliton 
part can be calculated explicitly. The crossing terms are close to the orthogonality condition (45) and 
hence are small. A combination of the quadratic terms in v,w gives the energy E (t) defined in (20), 
which is positive definite by Proposition 2. We thus end up with an estimate for ||w||l 2 + IMIff 1 if we 
can further show that H, U k , Q are almost conserved. 
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We first consider the angular momentum H k (t). Using the decomposition (41), we can show that 

Ilfc(i) - / p-% <d t ^d k <j>> da 

= / Pl\e < {g^dt^d^ > + P ; 1 q e (g £ ) M < 0^,0^ > da 

= [ p7W)°° < TPs + e ie (p e d e )- 1 w,d k ( ( f>s +e i@ (q e d e )- 1 v) > 

+ P ; 1 q e (gT < d^s + e i6 (q e d e )- 1 v),d k (<p s +e ie (q e d e y 1 v) > da. 

Recall that h G C 4 and h — m along the geodesic (t, uot). Since g e is close to h e in terms of the condition 
(55), we compare the above integral with that associated to the metric h. We hence can write Ilfc(i) as 
a sum of main terms plus an error 

n fe (i) =- < ip Sl d k <p s >dx - < ^ Sl d k {e m v) > dx - < e ie w,d k <j) S > dx - < e tS w, d k (e te v) > dx 
+ < (p e q € )- 1 (g e ) ol di(e i@ v),d k (e i@ v) > dx +Err(n k ), 

in which, by using Lemma 8, we can show that the error term Err(H k ) can be bounded as follows 

\Err(U k )\ < e 2 (l + |C| 2 )(1 + \\v\\ Hl ) + \\dg*\\ L ~(\\ W \\ L 4v\\ Hl + \\v\\ 2 H1 ) 
<(l + \tf)e 2 + e\\ W \\ 2 L2 +e\\v\\ 2 H1 . 

We show the crossing terms are vanishing due to the orthogonality condition. In fact, recall the definition 
of z in line (14). We find that 

dz dz 
D^ s = D z (p s ■ ^£ = -D z <p s ■ 7^ = -V^s. 

Since X(t) solves the modulation equations (51) and hence the orthogonality condition (45) holds, inte- 
gration by parts implies that 

< V>S, V x (e l6 v) > dx + < e ie w,V x <ps >dx= - < e i@ w, D^s >dx + < D s ips,e i@ v > dx = 0. 
We now compute the soliton part. Since dz = pdx, we can compute 

< ips, ^x4>S >dx =< e ie (ipujf u - p(u + u Q ) ■ V z f u ),e %e {-ipuf u {u + u Q ) + V z f u ■ ^) > dx 



= I P 2 u 2 fl{u + M ) + p(u + u ) ■ V z f u S7 z f u 
Jr 3 

= -pB(u + u ), 



dz 

TT dx 
ox 



where by Theorem 2, the ground state f u is spherical symmetric and B is given in line (57). We hence 
can write 

n fc (t) = P B(u k +u k )- < e m w,d k (e* e v) > dx + < (p e q € )- 1 (g*) ol di(e i@ v),d k (e ie v) > dx +Err(U k ), (71) 
where u = (u 1 ,u 2 ,u 3 ), u = (uI,Uq,Uq) and the error term Err(H k ) satisfies (70). 

We decompose the charge Q{t) in a similar way. Recall that (3 — pj 1 q e - We can show that 

Q(t)= / p: 1 q € <i(g e ) 00 Ws + e ie (p € d € )- 1 w),<j>s+e i@ (q € d e )- 1 v> 

JR3 

+ p- 1 q f _ < zig^Ok^s +e ie (q e d e )- 1 v)^ s +e l& (q e d e y 1 v > da (72) 
= - < iips, <f>s >dx - <iw,v > dx - < iips, e ie v > dx - < ie lB w, <f> s > dx 
+ < (p € q e )- 1 (g e ) ok id k (e i@ v),e i@ v > dx +Err{Q), 
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where the error term Err(Q) satisfies the estimate 

\Err(Q)\ < (1 + + £ |M| 2 L2 + e\\vf Hl . (73) 

Now observe that Dgifis — i4>s,Dgips = iips- The orthogonality condition (45) together with integration 
by parts implies that 

< iip s ,e te v > dx + < ie ie w,(t)s > dx =< D e ^s,e i@ v > dx - < Dg(j)s,e m w > dx = 0. 

For the soliton part, we can compute 

< itpsAs >dx= / -pufldx = -uj\\f u \\ L 2. 

Therefore, we can write 

Q(t) = w||/„||i 2 - < iw,v > dx + < (p e q e )- 1 (g e ) ok id k (e i6 v),e i@ v > dx +Err(Q). (74) 

The error term Err(Q) satisfies (73). 

Finally, we consider the main energy H(t), containing of quadratic terms and a higher order term 
corresponding to the nonlinearity. The quadratic part can be decomposed as follows 

^fe j. a j. ^ ^ aU a j ^ i™2 



< d k 4> > da - < d f (f>, d t 4> > da +m z <<j>,<j> > da 

'da 



-< {g e ) kl d k <f>,di(i> > da - < (g e ) 00 dt4>,d t <j) > da +m 2 < 4> s + e ie (g 6 d e ) _1 i;, 4> s + e l6) ( 9e d £ ) _1 w >c 



= I \V x <t>s\ 2 + \^s\ 2 +m 2 \<j> s \ 2 dx + m 2 < q; 2 v,v > dx + < q^(g*) kl d k (e* e v),d l (e* e v) > dx 
Jr 3 

+ 2 < V x (f)s,V x {e ie v) > dx + < w,w> dx +2 < e ie Wl ip s >dx +2m 2 < cj> s ,e ie v > dx +Err(Hq), 
where 

\Err(Hq)\ < (1 + |£|V + e\\w\\ 2 L1 + e\\vf m . (75) 
We expand the nonlinear term up to second order 

- / N{cj>)do = - [ N(f u + (qed^v) - N(f u ) - DNif^d^v - \{q e d c )- 2 vD 2 N{f ul )vda 

JR3 JR3 1 

-[ ^^ 1 + /£(9edO-V + ^/r 1 (« £ de)- 2 (H 2 + (p-l)K| 2 ) da 
Jr3 P + l ^ 

= - / -rrT/£ +1 + />i + J/rHkl 2 + (P - l)l«i| 2 ) + £rr(ffg), 
where by using Lemma 7, we can control the error term 

\Err(Hn)\ < \\vf Hl + \\v\\^ + (1 + |£| 2 )e 2 < \\vf Hl + (1 + |£| 2 )e 2 . (76) 
Group these together. We end up with 

n(t) = \ [ |V^s| 2 + Hs\ 2 + m 2 f 2 - ~^-Z +1 + m 2 q- 2 \v\ 2 + < q- 2 (g<) kl d k (e* e v),d l (e e v) > dx 

Z JR3 P + 1 

+ < V x (f> s ,V x {e i@ v) > dx +-<w,w > dx + < e ie w,^ s > dx +m 2 < <f> s ,e ie v > dx +Err(Hq) 
-<fZ,v 1>dx -\ f fP-^ + ip-^fP^vldx + ErriHn). 

£ ./»3 
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We can compute the soliton part 

1 JR 3 V + 1 



-5/ 

^ JR 3 



2pV/ 2 + (m 2 - c 2 )/ 2 + (± + ^ 2 )|V Z /| 2 - ^/r 1 ^ 



=P (^II/ W ||| 2 + ^I|V 2 / W ||| 2 ) =pB. 

Using the identity (16) and the orthogonality condition (45), we can eliminate the crossing terms 

< V x (f> s , V x {e ie v) > dx + < ^ s ,e ie w > dx + < m 2 (f> Sl e ie v > dx - < f*, Vl > dx 
= < fa, e m w > dx + < -A x <p s + m 2 4> s - e* e /P, e l& v > 
= < D x cf) S ■ V{\),c m w > - < D x iPs ■ V(X),e te v >= 0. 

Combining all these together, we can write 



H(t) = P B + l [ m\- 2 H 2 + < qj 2 {g*) kl d k {e i@ v),d l {e i@ v) > +\w\ 2 - 
/r 3 



2 



- (P ~ ^fZ^vldx + Err(Hn) + Err(Hq), (77) 
where the errors terms satisfy (75), (76). 

We proceed by arguing that the energy of the radiation term (v, w) is positive definite. Using the 
Fermi coordinate system, we have turned the stability of solitons on a slowly varying background into that 
on a small perturbation of Minkowski space. In Minkowski space, observe that the associated quantities 
"H, II fe , Q, as functionals of the solution of (7), are conserved and that the soliton (j) S is critical point 
of the Lagrange 

H - u k Il k - -Q 
P 

with Lagrange multipliers u k , — . To study the stability of stable solitons, we expand the above Lagrange 
around <f>s- It turns out that the soliton part is convex function of lo, u. The crossing term is vanishing 
since <j>s is critical point. We remark here that this is also given by the orthogonality condition (18), 
which has been shown above. The quadratic part gives the energy of the remainder <f> — <ps, which is 
proven to be positive definite under the orthogonality condition. One thus can obtain the orbital stability 
result of Theorem 3 in Minkowski space, see the work of D. Stuart [25]. We use a similar idea in our case 
by considering 

U(t)-(u k (Q)+u k ).U k (t)--^Q(t). 
Having the decomposition formulae (71), (74), (77), we can group the soliton part 

pB(l - (u(0) + u ) ■ (u + u )) - —\\f u \\l am . 

We denote the quadratic part 



E(t) = \ f 



™ 2 <Z7 2 M 2 + < qZ 2 (g*) kl d k (e ie v),d l (e i6 v) > +\w\ 2 - 



- (P - l)/<rS 2 + 2^ < iw, v > -2-^- < ( Pe q € )- 1 (g*) ol id l (e i@ v),e i@ v > 

+ 2 < e ie w,d k (e ie v) > (u k + u k ) - 2(u k (0)+u k ) < (p £ q t )- 1 (g^ al d l (e* e v),d k (e I0 v) > dx 
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as the energy of the radiation term (v,w). Since \u(0) + u \ < 1, p(0) > 1, \co\ < to, using Corollary 3 
and condition (55), we can conclude that 



E(t) >\( m 2 \v\ 2 + |V x (c^ e ^)| 2 + M 2 - fST>\ 2 " (P - + 2- < « > 

+ 2 <e ie W ,a fc (e 40 w) > ( U fe + ^)da;-Coe||i;||^ 1 
= ^ [\\w + p(u + u )V z v - ipuv\\ 2 L 2 {dz) + < V!,L + V! > dz + < v 2 ,L-v 2 > d ^j - C e||?;||^i 



(78) 



2p 



E Q (t) - C Q e\\v\\ 2 H1 , 



where the constant Co depends only on h, to. In the last line, we must view v, w as functions of z instead 
of x. The operators L + , L- arc defined in (19). We hence can write 



H(t) - (u k (0) + u k ) ■ U k (t) - -j--Q(t) =pB(l - («(0) + u ) ■ (u + «„)) - 

P(0) 

+ £(*) + £Vri + £rr 2 , 

where Srri consists of errors in H(t), II fe (t), Q(t) satisfying the estimate 



p(0) 



u\\L 2 



\Err-, I = 



Err(Hn) + Err(Hq) - (u k (0) + u^)Err(U k ) 



p(0) 



£rr(Q) 



^IHIin + (i + ieiV + e|Hilri + e|Hlii. 

Err 2 denotes the errors from the coefficients u + u , — , obeying the estimates 



|£>r 2 | 



-{u k - u k (0)) < e ie w,d k (e i@ v) > dx ~ ( 



~p ~ W) 



< IW, V >dx 



< 



u(t)-u(0)\\\w\\ L 4v\\ H1 . 



Now denote 



dH(t) = pB(l - (u(0) + uo) • (« + mo)) 



and 



Srr 3 =H(0)-(« fe (0) + ? 4)n fc (0) 



p(0) 
p(0) p(0) 



p(0) p(0) 



(Q(0)- W (0)||/J|| 2 (0)). 



We obtain 



B(t) + dW(t) =«(t) - w(o) - (u fc (o) + 4) • (n fe (t) - n fc (o)) 

(<9(*) - <9(0)) - Em - Err 2 + Err 3 . 



p(0) 



Since 



p(0)B(0)-| U (0) + no| 2 p(0)S(0) 



g(0) 
P(0) 



= 0, 



by Lemma 4, we can show that 



\Err 3 \ <e 2 . 



(79) 



(80) 



(81) 



(82) 



(83) 



Since (v,w) satisfies the orthogonality condition (45), Proposition 2 implies that Eo(i) is equivalent to 



the energy |H|^i 



when e is sufficiently small. Hence from (78), E(t) is equivalent to 



l_f/!(R3) 



L 2 (R 3 ) 



up to an error. The right hand side of (82) will be proven to be small in the next 



section. To obtain estimates for \\v\ 
fact, we can show 



2 2( R 3), it suffices to show that dH is positive definite. In 
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Lemma 9. Let X(t) G Ag . If Sq is sufficiently small, depending only A(O), then 

dH(t) > c(|w(t) - w(0)| 2 + \u(t) - u(0)\ 2 ) 

for some positive constant c depending only on A(0). 

Proof. As functions of u, we can compute 

D u (p(l - (u(0) + n ) • (u + uo)))(0) =D P (0)p(0)- 2 - p(0)(«(0) + u ) = 0, 
D£(p(l - («(0) + m ) • (« + uo)))(0) =D 2 p(0)(l - |u(0) + u | 2 ) - 2£p(0) • (u(0) + «o) 

=p(0)7 + p(0) 3 (u(0) + u Q ) • (u(0) + uo), 

which implies that the Hessian of p(l — u(0) • u) with respect to u is positive definite. Hence, if S n is 
sufficiently small, depending only on A(0), then 

P (l - (u(0) + uo) • (u + mo)) > piO)- 1 + d\u(t) - u(0)\ 2 , \u(t) - u(0)| < (5 

for some constant ci depending only on u(0) + wo- 

For the other part on the right hand side of (81), we rely on the properties of the ground state f u . 
Recall the definition (57) of B(t) and the scaling property (11) as well as the energy identities, we can 
show that 

D w (B(t) - o; 2 ||/ w ||| 2 - B(0) + w(0)w||/ u ||i a (0)) (0) = 0, 
Dl (B(t) - Lo 2 \\.U\ 2 L2 - B(0)+l,(0)lu\\.UU0)) (0) 

= (^Y-(O) 2 - m 2 ) (m 2 c(0) 2 )«||.f \\ 2 L2 > 0. 

The positivity is due to the assumption A(0) G A^. Hence if So is sufficiently small, we can conclude that 

B(t) - lo 2 \\.U\ 2 L 2 - B(0) + lo(0)lu\\.U\ 2 l2 (0) > c 2 \u{t) - Lu(0)\ 2 ycu(t) G (u - 5 q ,lo q + S ) 

for some constant c 2 depending only on A(0). 
Therefore, we have shown that 

dH(t) > piOy'Bit) + Cl \u(t) - u(0)\ 2 B(t) - p(0)- V 2 |l/J! 2 + B(0) - c(0) W ||/ w || 2 2 (0)) 
> ci|«(t) - u{i))\ 2 B{t) + c 2 p(0)- 1 \uj(t) - u(0)\ 2 > c(|w(t) - lo(0)\ 2 + \u(t) - u(0)\ 2 ) 

for some constant c depending only on A(0) if 5q is sufficiently small. □ 

We now choose 8q such that Lemma 9 holds and the set A«5 defined in line (52) is nonempty. By 
Proposition 2, we have shown the left hand side of (82) is bounded below as follows 

E(t) +d%> c(\\v\\ H i(t) + \\w\\ L 2(t) + \u(t) - co(0)\ 2 + \u(t) - u(0)\ 2 ) 

for some constant c depending only on A(0) if e is sufficiently small. 

6.4.4 Energy Estimates 

We estimate the right hand side of (82) in this section by using the integral formulae (64), (66), (69). 
Denote the soliton part 

/ / T» v [<l>s]rfzdxdt, a G {0, 1, 2, 3}, 
Jo Jm 3 

where d t (j)s should be replaced with ips- 
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Proposition 6. We have 



\H(t) n(0) + H° s (t)\ <e 2 + ef \u\ 2 + \\vf H1 + \\w\\Us, 

Jo 

\u k (t) n fe (o) - u k s (t)\ <£ 2 + tf M 2 + \\v \\ 2 H , + \\M\l*ds, 

Jo 



IQW-Q(0)l<e 2 + e / 

Jo 



M 2 + IMI#i + IMIi 2 ^ s 



if e is sufficiently small. 

Proof. We decompose the energy momentum tensor T^tfi] as follows 

Tfiui^] = Tfivi^s] + T^[(j) S , 4> - 4>s] + QT^ v [(f> - 4>s], 

where 

T^sA- <l>s\ = < d^s,d v (^)- (p s ) > + < d v (j)s,d^((j)- 4> s ) > 

- 9%{< d^s, - <p s ) > + < V'(<fe), <P~4>s >), 

QT^ u [(p—(j)s\ is at least quadratic in v, w. Again, dt4>s should be replaced with ipg in the above expressions 
as well as in the following argument of the proof. Recall the deformation tensor (n da )^ v = — \d a (g € Y v . 
Using Lemma 7 to estimate the nonlinear terms in QT^cj) — <fis}, we can show that 

f j T^-<j> s }(ir d «r»dvol <e [ M^ + Mlids, Va e {0, 1,2,3}. 
Jo Jm 3 Jo 

For the linear term T M „[</>s, </> — <j) S ], we apply Lemma 8 to get 

f ( T^[<t> s ,<t>- <j>s] (7r a «)^rfvol <e 2 f\\\v\\ H i +\\w\\ L *)(l + \ti(s)\)ds, Vae{0,l,2,3}. 
Jo Jr 3 Jo 

For the soliton part, first we can drop the volume factor d 2 as 



o Jm 3 



T^s^ndl-^dxdt 



< , 3 



f\+\as)\ 2 ds. 

Jo 



The above estimates rely on the unknown upper bound of Although it is bounded by Ci as a 

boostrap assumption, we do not want C*2 to appear in the following estimates. Otherwise, we need 
smallness on t in order to close our boostrap argument. We will instead use GronwalPs inequality to 
control |£|. First, we use the relation £ = u + f) together with Corollary 4 to control £ in terms of u, v, 
w. We can control f] by Corollary 4 as follows 

\v\ 2 <\M 2 <(C2e) 4 + \\vf H1 <e 3 + \\vf H1 . 

If C|e < 1, then we can show that 

\m? < \m\ 2 + t f \u\ 2 + \i\ 2 ds < 1 + e- 1 f \u\ 2 + \\vf H ,ds, Vt < Tje. (84) 
Jo Jo 

Plug this into above estimates. We conclude that 

\U(t) - H(0) + U° s (t)\ < 6 3 f 1 + \tfds + e f \\vf Hl + \\ W f L2 ds 

Jo Jo 

<e 2 +e f\u\ 2 + \\vf Hl + \\w\\l>ds. 
Jo 
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Similarly, we can show that 

U k (t) - n fe (0) - H k s (t) - f ( {rn Q y v d v pT™°[<t> s ]dxds 

Jo Jts.3 



Q(t) - Q(0) 



o Jm 3 



< idnf3(ma) ,lv d v <f>s, <t>s > dxds 



<e 2 +e 



<e 2 +e 



f 

Jo 

f 

Jo 



\u\ 2 + \\v\\ 2 H1 + \\ W \\ 2 L2 ds, 



+ w 



L^ds, 



where we recall that mo is the Minkowski metric, (i = pj 1 q e and we denote T™°[(j)s] as the energy 
momentum tensor associated to the Minkowski metric mo. We thus can conclude the proposition if we 
can control the two integrals 



(moY v d v f3T%°[<i>s}dxds 



< idf,/3(ma)^d u <j>s,(ps > dxds. 



The only smallness in the above integrals is contributed by <9/3, which has size e. To prove that they are 
higher order error terms, we have to exploit the properties of the solitons 0s- The key observation is that 
the solitons travel in the direction dt + u^dk- More precisely, we compute the second integral 

< id^(m r v d v <f> s , <f>s >= -d t P < *Vs, <Ps > +dkP < id k <Ps, 4>s >= P^f 2 (z)(d t + (u + u )V)/3. 

Observing that z — A u (x — £ — u t), f u (z) travels along the geodesic (t,uot), thus integration by parts 
may allow us to gain extra smallness. For any function F(z, X(t)) independent of O, we have the identity 

{dt + (u + u )V x )po ■ F = d t (p F) - p D x F ■ (V(X) +j) + (u + u )V x p ■ F 

= d t (J3 F) - fioDxFj - p (u + u )D i F + (u + u ) V x /3 ■ F (85) 
= d t (0oF) - (3 D x Fj +(u + u )D x {p F), 

where we have used D^F = —D X F as z — A u (x — £ — u$t). This key observation (85) is of particular 
importance in this paper. It allows us to prove Theorem 4 under the sharp condition q > 1. Now let 
/3q = f3 — 1, F = poofu(z). The inequality (85) yields the estimate 



/'/ 

Jo Jr 3 



< idu,f3{m a y v d v (f>s,(l>s > dxds 



<e 2 (l + |£(t)| 2 ) + | 7 |e 2 (l + |?( S )| 2 )ci S 



\ u \ + \\ v \\w ds ) + C 2 e 



<e 2 + e 



7 

Jo 

/V + IMI^ + Hli^ 

Jo 



2,2 



Jo 



2 £ ) 2 + IMIffi^S 



for all t < T/e by Corollary 4 and the boostrap assumption (53), C|e < 1. We thus conclude the first 
and the third inequality of this proposition. 

To show the second inequality of this proposition, it suffices to estimate the first integral above. 
Similarly, we can compute 

= {moT v d v ((/3 - 1)T^[4> S ]) - (/? - l){m Y»d v T™°[cl>s] 
= {m y v d„ ((/? - l)T™°[<i>s]) - 08- 1) < -d t ij s + A x cf> s - m 2 s + ^sK^sMs 
= (tnorflv (os - i)T^y> s ]) + 08 - 1) < D X ^ S ■ 7,a fc 0s 



> 



> 



by using the identity (16) and the relation A = V{\) + 7. Hence according to Lemma 8, we can show 
that 



{m Q y v d v pT™°[<]>s]dxds 



Jm 3 



We hence have proven the proposition. 



< 



e 2 (l + 



<e 2 +e 



Jo 

f\u\ 2 + \\v\\ 2 m + \\w\\l>ds. 
Jo 



□ 
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This proposition allows us to control the right hand side of (82). 
Corollary 5. We have 



H(t) u(0) (u k (o) + u k Q ) ■ (n fe (t) - n fe (o)) - -^m) Qm 



<e 2 



<-- ! I \A 2 + M\m +\\wf L *ds. 



Proof. Notice that T^[(f>s] is a function of (z, X(t)) and is independent of 9. By applying (85), we can 
show that 



THMd* + (« + u Q )V x ){g e )^dxds 



10 JR 3 

Since tt^° = \d a {g e )^ v , we have 

\H° s (t) + (u(0)+uo) k H k s (t)\< 



+ 



<e 2 +e [ \u\ 2 + \\v\\ 2 Hl + \\w\\ 2 L2 ds. 
Jo 



T^[<f> s ]{dt + (u + u )V x )(g e )^dxds 

( f T»»[<p s ](u(s) - u({))) k d k (g')^dxds 
Jo Jr 3 



<e 2 + e \u\ 2 + \\v\\ 2 H1 + \\w\\ 2 L2 d S + e 2 (1 + + e)ds 



<^ 2 + e /V| 2 + M| 2 ffl + IIHI| 2 ^ 
Jo 



where we have used \u(s) — u(0)\ < \u(s) \ + e by Lemma 4. Then the corollary follows from Proposition 
6. □ 

6.4.5 Proof of Proposition 4 and Theorem 5 

We are now able to improve the boostrap assumptions and to conclude Proposition 4 and Theorem 5. 
Denote 

£{t) := M 2 + Mt)- w (o)| 2 + |HI 2 L2 + IMI 2 ff i- 

According to Proposition 2 and Lemma 9, we have 

E{t)+dH(t) > c{l-C e)S(t) 

for some positive constants c, Co independent of e, Ci- Hence by (82), (79), (80), (83) and Corollary 5, 
for sufficiently small e, we can show that 



m < e 2 (l + ICI 2 ) + e f £(s)ds + £{t)i + e£{t) 
Jo 

<e 2 + e j £{s)ds + £{t)i +e£(t), 
Jo 



where we have used inequality (84) to control |£| 2 . Since the implicit constant is independent of e and 
initially £(0) < e 2 , using Gronwall's inequality, we have 

£(*) < e 2 , Vt < T/e. 

Let C3 be the universal implicit constant appeared before. By our notation, C3 depends on h, m, A , T 
and is independent of e, C2. In particular, we have 

£(t) = \u(t)\ 2 + \oj(t) - oj(0)\ 2 + \\v\\ 2 m {t) + \\w\\ 2 L2 (t) <C 3 e 2 , Vt < T/e. 
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By Corollary 4, this implies that 



\m\< / \u\ + n\ds 




<C 3 e + C 3 e + C 2 e 2 + \\v\\ 2 H1 ds 



Jo 



<C 3 e + cjT + C 3 T + C'leT, 



where we let e to be small such that Cle < 1. Take 



C 2 =C 3 + cjT + C 3 T + 2C 2 T. 



If 



e^min^- 4 ,-^- 1 ^,!}, 



then for all t < T/e, we have 



|£| < C e + c|t + C T + C 2 eT < C 2 . 



This improves the boostrap assumptions (54), (53). Therefore we can conclude that Proposition 4 follows 
from Corollary 4. 

For Theorem 4, notice that there is a unique short time solution <p(t, x) on [0, t*) x R 3 . Then Lemma 6 
implies that the modulation equations (51) admits a local solution X(t), t G [0,t**] C [0,t*). Proposition 
4 then shows that <j>(t,x), X(t) satisfy (39), (40) with a constant C independent of e. Since Proposition 
4 holds for any t** < T/e, we can conclude that the solution <p(t,x), A(i) can be uniquely extended to 
[0,T/e] x M 3 and satisfy estimates (39), (40). 

6.5 Estimates of Higher Sobolev Norms 

In the previous section, we have proven the orbital stability of stable solitons of equation (36) on a slowly 
varying background in the energy space H 1 x L 2 . To solve the full Einstein equations (1) and to obtain 
a C 1 spacetimc ([0,t*] x R 3 ,g), we need higher Sobolev estimates for the matter field <p. In this section, 
we prove Proposition 3 and Corollary 1. 

Since we already have a solution <p(t, x) and a curve X(t) satisfying estimates (39), (40), we use energy 
estimates to obtain the higher Sobolev estimates by considering the equation of the remainder v. However, 
to simply the argument in the sequel and to avoid taking fourth order derivative of "f(t), we choose a 
modified curve X(t) € A sta b, defined as the integral curve of V(X), that is 



Then we claim that Proposition 3 is reduced to the following estimates. 

Proposition 7. Assume the initial data </> <= H 3 , e H 2 . Let e\ be defined in Proposition 3. Then 



d t x = v(x), x(o) = (w ,e ,o,u ). 



Using this modified curve X(t), we decompose the solution (j> as follows 



<t>(t,x) = 4> s (A; x) +c 



(86) 



]T \\d s v(t,-)\\ L 2 {R3) <max{ e , ei }, Vt<T/e. 

\s\<3 
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In fact, if Proposition 7 holds, observing that 

|A-A|< r|7l^+|A(0)-A(0)|<6 2 t + 
Jo 



T 

e<£, Vt<-, 



<9 t 6(A) = ~ + Powo«o(mo + u), V X @(A) = -p u} u , 



then we have 



||0 - 4> s (x; X)\\ H s = Us(x; A) + e ie ^v - 4> s (x; X)\\ H s 
< \\(t)-\(t)\ + \\v\\ H 3 <max{e,ei}. 

This partially explains Proposition 3 as ||<9 t s (0 — 0s(x; A))||^2 requires taking higher order derivatives of 
the modulation curve X(t). We will estimate the higher derivatives of X(t) in next section. 
Although we have modified the curve X(t), the remainder v is still small in H 1 . 

Lemma 10. We have 

Nlifi + \\dtvh' < e, Vt€[0,T/e]. 
Proof. Since \\v\\ H i + ||w>||l2 < e, by using the decomposition (41), we can show that 

\\e mCx) i\\m = ||^(A;x) + e U) M{q e d e )- 1 v - Mk*)\\m < C 
|| e z©CA )6 , t ~|| i2 = ||a t (0_0 s (A;x))- J 9 t e(A)e te ^5|| L 2 

= ||D A 0s(A; x) • V(X) + e^&d^w - D X ^ S (X; x) ■ V(X) - *5 t e(A)e* e (%|| L2 < e. 

□ 

Remark 8. TTie reason that we consider the modified curve X is that we must avoid taking the fourth 
derivative ofX. Otherwise, we have to take third derivative of the nonlinearity which is impossible 

since p is assumed to be less than 3. Using the modified curve guarantees that when we differentiate the 
equation of the radiation term v twice, we only need third derivative of X. We remark here that X still 
depends on X. 

6.5.1 Estimates for Higher Derivatives of X(t) 

Since A(i) satisfies the modulation equations (51), by differentiating the equations, we are able to derive 
estimates for derivatives of X(t). It turns out that we first have to estimate the derivatives of the 
nonlinearity Af(X). 

Lemma 11. Let AT (X) be defined in line (49). Assumep>2. Then for any vector field Y on [0,T/e]xK 3 , 
we have 

WAji^M + K-^dy^i + iyvi + H), 

\Y 2 ^W\ <(1 + NU~ + \Y ln/ w |)(|Y/ w | 2 + \Yv\ 2 + \v\ 2 ) + (\v\ + M^XI^/d + \Y 2 v\) 
+ \Y 2 d £ \(\v\ 2 + \v\P). 

Proof. It follows by direct calculations and the properties of f u summarized in Theorem 2. □ 

Using this lemma together with the modulation equations (51), we are able to estimate the higher 
order derivatives of the modulation curve X(t). 

Proposition 8. Let j(t) be defined in line (17). Assume -f(t) satisfies the modulation equations (51). 
Then we have 

l7l<e 2 , 

\dh\<e 2 (l + \\v\\ L ~) + 6110^X^11^. 
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Proof. Differentiate the modulation equations (51), we get the ODE for 7 

(D + A + D 2 )j + d t (D + £>! + L> 2 ) 7 = 5 t F(i; A(t)). 

Since |A(t)| < 1, \dg e \ < 1, we can show that 

|5 tJ D| + |9 tJ Di|<|A| + |9 5 e |<l 

according to the definition of D, D\ given in line (50). For d t D 2 , we have to use the equation of <p. In 
fact, by (36), (47), we can show that 

\d t D 2 \ < |A| • \\w\\ L 2 + I < D\^s, d t (e ie w) > dx \ + |A| • \\v\\ L 2 + \ < D 2 x <f> s , d t (a (<j> - fa)) > dx \ 

< 1 + I < D\dp s , {<t>tt ~ d t ^s)ai > dx \ + \< Dfys, a (d t <t> - D X ^ S (V(X) + 7)) > dx \ 

< 1 + I < Dl<jy s ,H(t, x)+A^- b(m 2 - V) >dx I + I < Dj^s, ao(<k ~fa~ Dxfai) >dx 

< 1 + I < d k {a k ^D\<j)s), > dx \ + \< V x Dl4> s , V x </> > dx \ 

< 1 + \\v\\ H i + \\w\\ L 2 < 1. 

We hence have shown 

|$(£> + £>i + £> 2 )| < 1. 
Then Lemma 6 implies that estimates \ j\ < e 2 follow if we can show that 

|d t F(t;A(t))|< e 2 . (87) 

By the definition (48) of F(t; A(t)), it suffices to estimate 

< bDxfa, e l@ d t Af{\) >d x , < diDxfa,fa >d x , < d'oDxfa, fa- fa >dx ■ 

All the other terms can be estimated similarly to F(t; A(t)) in Lemma 5. We first consider the nonlinear 
term, which is in fact the main term in dtF(t; X(t)) as other terms are errors from the slowly varying 
metric g e . The key observation is that for vector field X = dt + uq\7 x and any C 1 function F\, we have 

|XJi(/ w (z))| <\(8 t + u V x )Uz)\ = iD^f^io + \7J u z + u \7J u A u \ 

< \M + \^zU{A u u + A u (-i - u ))\ (88) 

< e 2 + |«| < e 

by using Theorem 5 as well as the relation £ = u + f\. Hence by Lemma 11 and Lemma 7, we can show 

I < bD x <t> s ,d t N{\) > dx \<\< bD x <t> s ,XAf{\) > dx I + I < bD x <t> s ,u V x N{\) > dx \ 

< \\D x <t> s (\v\ + Iv^QXfvl + \Xv\ + \v\)\\ L i + I < u \7 x (bD x <t> s ),Af(\) > dx \ 

< (\\D x <j> s ■ Xv\\ L 2 + e)(\\v\\ L 2 + Wv^- 1 ) + \\v\\ 2 Hl + \\vf H1 

<e 2 + e\\D x 4>sd t v\\ L 2 < e 2 + \\D x <j> s d t (e ie v)\\ L 2 + \\D x <j> s ev\\ L 2 (89) 
<e 2 + e\\w\\ L 2 <e 2 . 

We must remark here that using the decomposition (41) we have 

\d t v\ < \w\ + \Qv\. 

Although 9 depends on z, \QD X (j>s\ decays exponentially in z by Theorem 2. 

For the second term < a\D x 4>s,fa > dx , we use equations (36), (47) and then use integration by 
parts. We can bound 

I < cuD^s, fa > dx \ = \<d t In ai D x ^ s , a^d^ + 6"0„0 + - b{m 2 - \<f>\ p - V) > dx \ 
<e 2 + \<d k {a k > i d t lna 1 D x <j )s ),d^> dx \ + \<\7 x (d t \na 1 D x fa),\7 

<P > dx I 

< e 2 + \\d 2 g^D x fa\ L 2 <e 2 + ||c>V - h c )\\ L , + \\d 2 h^D x fa\ L 2 < e 2 
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by using the assumption (37) and the fact that h e (t,x) = h(et, ex) <E C 2 . 

The third term < d'oD x ips, <j> ~ 4>s >dx can be estimated similarly. We can show that 

| < doD x il>s,<j>- <fe >dx I < e 2 + \\d 2 g e D x <t> s \\ L 2 < e 2 - 
Therefore, we have shown (87). Hence \^\ < e 2 . 

Having proven \^\ < e 2 , to estimate of the third order derivative of 7, differentiate the modulation 
equations (51) twice 

(£> + £>i + D 2 )dh + 2d t (D + D,+ D 2 )j + d tt {D + D 1 + D 2 )j = d tt F(t; A(t)). 

Similarly to estimating dtF(t: X(t)) carried out above, we can show that 

\d tt (D + D 1+ D 2 )\< |A| + |A| + \\d 2 (g e - h*)\\ L2 < 1. 

It hence suffices to estimate duF(t; X(t)). The strategy is similar to that of dtF(t; A(t)). The main term 
is that with derivative hitting on the nonlinearity < e~ ie bD X (fis , d t tN(X) >, which by Lemma 11 and by 
inequalities (88), (89), can be estimated as follows 

\<e- m bD X (b s ,d tt Af(X) > dx \<\< e-^bD^s^duAfW > dx \ 

< I < e- i@ bD x <f> s ,X 2 Af(\) > dx \ + \< e- ie bD x <j> s , (-(u V x f + 2u Q \7 x X)Af(X) > dx \ 

< e 2 (l + ||w||ioo) + e\\D x cj> s X 2 v\\ L 2 + \\D x <p s d 2 g^\v\ 2 + \v\*)\\ L i 

+ I < (u V x ) 2 (e- i@ bD x <j> s ),M(X) > dx \ + \< u Q \7 x (c- ie bD x cf>s),XAf(\) > dx \ 

< 6 2 (1 + \\v\\ L oo) + e\\D x cj>sX 2 v\\ L 2 + \\d 2 (gt - h*)\\ L 4v\\ 2 H1 
<e 2 (l + \\v\\ L -) + e\\D x <p s X 2 v\\ L 2. 

For those terms when the derivative hits on dt<p, we rely on the equation (36) together with the identity 
(47) and then use integration by parts to pass the derivative to the metric g e . We hence can show that 

\dttF(t; X(t))\ < e 2 (l + \\v\\ L ~) + e\\D x s X 2 v\\ L , + \\\d 3 (g<)\(\D x <b s \ + \D x 1> s \)\\v 
<e 2 (l + \\v\\ L ~) + e\\D x s X 2 v\\ L2 . 

Then Lemma 6 yields the estimate 

\d?l\ < e 2 (l + Nk~) + e\\D x s X 2 v\\ L 2. 

□ 

6.5.2 Linearized Equation for v and Energy Estimates 

Using the modified curve X(t) and the corresponding decomposition (86), we can find the equation for v 

L e v + Af(~X) + F = 0, (90) 

where 

F = C - ie ^{U gC s (x-x) - m 2 (f> s + |<fe| p_ Vs) + «<r© ' v. 
For any complex function v(t, x) — Vi(t,x) + iv 2 (t,x), the linear operator L e is defined as follows 

L e v = U g ,v + A(i)v + 2id^e ■ d»v - m 2 v + fZ~\z)v + (p - l)/^" 1 ^)^ (91) 

with 

= -{gT v dMl)dMl) = e- m a ge e ie - in fl .0. 
We have the following energy estimates for the linear operator L e . 
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Lemma 12. For all t < T/e, we have 

II^IM*) 2 < ll^(0,x)||| 2 + 1^(0,^)111, + e" 1 /' \\L e v\\Us)ds + s 

Jo o< 



sup \\v\\ L2 . 

<s<t 



Proof. Recall the energy momentum tensor T^u] for the operator D s 



Tn V [v] =< d^v,d u v > --g% < <9 7 v,<9 7 « > . 



For any vector field Y, we have the identity 

r>"(r„„Hy) = t<">]<,+ < n^Ytf) > . 

Take Y = X = d t + u^dk- Then integrate on the region [0, t] x I 3 . We obtain 

/ f^[v]XVda(t)= I f flv [v}X»n v da(0) + I I f^[v]n^+ < U g ,v,X{v) > dvol. 
jr 3 Jr 3 Jo Jr 3 

By replacing O g ,v with L e v, we must estimate the other terms respectively. First, consider the term 
< 2id^Q ■ d^v,X{v) >da- We use integration by parts. We can write 

< 2id fi e ■ d^v, X(v) >d 2 = d^Q < 2id fi v, Xv >= d\d»Q {d^ < iv, Xv > -X < iv, d^v >) 

= cV (d 2 d"Q < iv, Xv>)-X (d 2 d»Q < iv, d^v >) - d^d^Q) < iv, Xv > +X{d 2 d»Q) < iv, d„v > . 



Recall the definition of A. We can compute 



d t <d = h p oj (\u \ z + u u) = p uj + p uj u u - 

P P Po 

V 2 = -PqUqUq. 



Therefore 



< 2id u G ■ d»v,Xv > dvol 



JR 3 



< 2id fl Q ■ d»v,Xv > d 2 dxds 



< 



< 



I <iv,d°eXv-d^ed„v>da +11 \d(d 2 d^Q)vdv\dvol 
Jr 3 o Jo Jr. 3 

v\\ L 2(t)\\dv\\ L 2(t) + ||«|| ia (0)||ft;|| L2 (0) + e I \\v\\ 2 L2 + \\dv\\ 2 L2 ds. 

Jo 



For the other terms, notice that for any real function F\, we have 

2 I I < F lV ,X(v) > dvol= I I F!X(\v\ 2 )(%dxdt 
Jo Jr 3 Jo Jr 3 



= I F^vfda -II X{F l d 2 e )\v\ 2 dxdt. 
Jr 3 Jo Jr 3 



Jo 

When Fx = fP' 1 , since |XFi| < e by (88), we have 

When Fi = ^(7), since \d®\ < 1, \d 2 0\ < |<y| < e 2 , we have 

WXiA^yh, = \\X{{fy»dp®d v ®dl)\v\ 2 \\ L , < e\\v 2 \\ L1 < e\\v\\ 2 L2 . 
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Combine all these together. We have shown that 

/ / <A( 1 )v-m 2 v + fP; 1 (z)v + (p-l)fP; 1 (z)v 1 ,X(v)>dvo\ 
Jo Jm 3 

< \\v\\ L2 (t)\\v\\ m (t) + |MM0)HM0) + e f \\v\\ 2 m ds. 

Jo 

Now, since the unit normal vector field (— (h e )°°)~^ (/i e )° M 9 M to the hypersurface R 3 as well as the vector 
field X — dt + Ugdk are timelike with respect to the metric h e , we can conclude that the unit normal 
vector field n = (— (g e ) m )~ 2 (g e ) 0fl d^ together with the vector field X — dt + u^dk are also timelike for 
the metric g e if e is sufficiently small. Therefore, there exists a positive constant c, depending only h,u 0} 
such that 

T^[v}X»n v >c\dv\ 2 . 
Since \n*\ < e, the energy identity then implies that 

\\dv\\Ut) < ||^(0,x)||| 2 + \\v(p,x)\\h + e f \\dv\\Us)ds + e" 1 f \\L t v\\l^s 

Jo Jo 

+ \Hh+e f ' \\v\\ 2 H1 (s)ds, Vt<T/e. 
Jo 

The Lemma then follows by using Gronwall's inequality. □ 
6.5.3 H 2 Estimates 

Commuting the equation (90) with the vector field X — dt + u^dk, we obtain 

L e Xv + [X, L c ]v + XAf(v) + XF = 0. (92) 

Since we have computed 

|36|<1, |d 2 9| < |7| < 6 2 , \8Am<e, \Xf u \<e, 
using Lemma 10, we can estimate the commutator 

\\[X,L e ]v\\ L , < || [X, U gC ]v\\ L , + \\X8»0 ■ d^\\ L 2 + \\Xf£\v + (p- + \\XA( 7 )v\\ L 2 

< e\\d 2 v\\ L 2 + \\d 2 (g e )dv\\ L 2 + e(\\dv\\ L 2 + \\v\\v) 

<e\\df i ^ L , + \\af i {g e -h e )\\ L e\\dv\\ La + ^ 

<e\\d 2 v\\ L 2 + \\d 2 r\\ H i \\dv\\ H i+e 2 

< e\\d 2 i\\ L 2 + e 2 \\dv\\ m +e 2 < e\\d 2 i,\\ L 2 + e 2 . 

For the nonlinearity XJ\f(\), Lemma 11 yields the estimates 

\\xMCx)\\^<\\(\v\ + \ir 1 )(\x.U + \xi\)\\ L 2 

<\\\v\ + \v\P-^\\ L3 (e+\\Xv\\ L e) 

< (\M& + WvWwXe + \\d 2 v\\ L z) <e 2 + e\\d 2 i,\\ L 2. 

As for XF, first using the identity (16) we have 

U4> s {x; A) - m 2 <f> s (x; A) + |0s| p_1 0s = V(\)D 2 ^ S V{\) - V{\)D 2 <p s V{\) - D x <j> s ■ d t V(X). 
By Theorem 5, we can show that 

\V(X) - V(X)\ = 1(0, -, u + u , 0) - (0, —,^o,0)| < e. 

P Po 
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The key observation is that for all k < 4, we have 

\\XV k U\ L 2 <e, 

which can be proven similarly to (88). In particular, we have 

\\X( C -^D 2 ^ s )\\ L ,<e. 
Therefore, by Proposition 8, we can estimate XF as follows 

\\XF\\ L 2 < \\X (v(\)c- ie Di0 s V(\) - V{\) C - ie Dl^ s V(\) - e~ iS D x <!>s ■ d t V(\)) \\ L , 

+ \\X(c-^(a gC - a)d> s )\\ L2 + ||X(D s .0 • «)|| La 

< |i| + |7l +e\\X(c-^Dlts)\\Li + \\d 2 (9 e )(\V 2 .U + NU+f^Wv 

+ ||(.g e -m )(|V 3 .U + |V 2 .U| + \VU+f U0 )\\ L 2 + \\d\g*)v\\ L 2 + \\d(g e )dv\\v 

<t 2 + \\d 2 r\\m\\v\\Hi+4dnLi<t 2 - 
The energy estimate Lemma 12 then implies that 

||3^|| 2 L2 < max{ £ 2 , £ 2 } + f \\[L e ,X}v\\ 2 L2 + \\XM(\)\\h + \\XF\\l*da (93) 

Jo 

<max{e 2 ,e 2 } + e / \\d 2 v\\ 2 L2 ds 
Jo 

for all t < T/e, where e\ is the size of the initial data given in Proposition 3. 

To derive the full estimates for ||i9 2 {;||^2, merely commuting the equation with the vector field X is 
not sufficient. The key point to obtain estimates (93) is that the soliton <j>s travels along the timelikc 
geodesic (t, uot) or quantitatively the vector field X = dt + uoV x acting on f u leads to the estimates 
\Xf u \ < e. For general vector field, estimates (93) may not hold. To retrieve the full estimates ||9 2 i;|| i 2, 
we rely on the following elliptic estimates. 

Lemma 13. Let A lJ (x) G C a (R 3 ) for some positive constant < a < 1. Assume that A 1 ? is uniformly 
elliptic. That is 3K such that 

K-^yl 2 < A^(x)y iyj < K\y\ 2 , Vx,y e R 3 . 
Then there exists a constant C such that 

M m <C\\A\\ 2 Ca (\\A^d lJ cj ) \\ L , + U\\ L2 ) 
for any (j) e H 2 (R 3 ). Here \\A\\ C <* = sup p ij '||c«. 

Proof. Let x be a cut-off function supported in the ball B2 with radius 2 and equal to 1 in the unit ball 
B\ . Then the elliptic estimates show that 

U\\ H 2 {Bl) < q|A|| ca (p^( x <£)|| L2 + H\\ L 2 { b 2) ) < c\\A\\ ca (\\A^d ij( t>\\ L 2 {B2) + U\\ HH b 2) ) 

for some constant C independent of <f), \- The above estimate holds for any ball B 2 - We cover the whole 
space M 3 with radius 1 balls such that every point is covered for at most 10 times. Add all the estimates, 
we conclude that 

H\\h> <C\\A\\ Ca (\\A^d ij( f>\\ L2 + U\\ Hl ). 
Interpolating between H 2 and L 2 , we have 

U\\ H1 < cwnUmh < Ic-'WAWclUWH* + 2C\\a\\ c -M\\v- 
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Plug this into the above inequality. We get 

UWm <q|A|| 2 c 4||^0|| i2 + ||0|| i2 )- 

□ 

Using this lemma and estimates (93), we are able obtain the H 2 estimates of v. First, write the wave 
operator D ge as follows 

= i(9 e ) kl + (g e )°°44 ~ 2{g f -f k u l )d kl i + ((g e f°(d t - u V) + 2{g*f k d k )Xi + d^d^MDdvV. 
Since X — d t + u k dk is timelike with respect to the metric g e , we conclude that the 3x3 matrix 

(g*) kl + (g*) 00 u k u l - 2(f) ok u 
is positive definite. That is there exists a constant K, depending only on h, uo such that 

K- l \y\ 2 < Y, ykyi((9 e ) kl + (9 e ru k u l -2(g-ru l )<K\y\ 2 , Vt < T/e, x,yeR 3 . 

l<fe,i<3 

Hence by Lemma 13 and equation (90), we can show that 

\\v\\ H 2 <e+\\dXv\\ L2 . 

Then inequality (93) implies that 

||<9X5|| 2 2 < max{e 2 ,e 2 } + e /" e 2 + \\dXv\\ 2 L2 ds < max{e 2 , e 2 } + e f \\dXv\\ 2 2 ds, Vt < T/e. 

Jo Jo 

By using Gronwall's inequality, we conclude that 

||9 2 5||| 2 <max{6 2 ,6 2 } + ||M5||| 2 <max{e 2 ,6 2 }, Vt < T/e. (94) 

6.5.4 H 3 Estimates 

Having the H 2 estimates, we first can improve the estimate of df-f obtained in Proposition 8. According 
to the decomposition (41) corresponding to the curve X(t), we can show that 

Ml- < INI* 2 < ||aoc- je(A) (^s(A; a; )-^(A;x)-e je W^)|| ff2 
<\\(t)-\(t)\ + \\i\\ H 2 + \\d 2 a i\\ L2 } 
< max{e,ei} + \\v\\m\\d 2 a \\ H i < max{e,ei}. 

We also need to estimate X 2 v, which, as having pointed out previously, does not follow directly from the 
estimates of v. However, notice that 

\\D^ s X 2 {a e-^)\\ L . < 1 + ||X 2 a || L2 < 1, 

||0s(A;a;) - 0s(A;x) - e^v^- < ||5|| LOO + |A(t) - \(t)\ < max{e,ei}. 
By the decompositions (41), (86), we can show that 

WDtfsXMv < \\D^ S X 2 (aoe-™(M**)-MKx)-e™ W{X) v)) IU 2 
< max{e, ej + |A(t) - A(t)| + |-y| + ftl < max{e, Cl }. 
Thus Proposition 8 implies that 

\dh\ <e 2 (l + \\v\\ L ^)+e\\D x ^ s X 2 v\\ L 2 <emax{e, ei }. (95) 
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We proceed to estimate the H 3 norm of v. Commute the equation (92) with the vector field X = 
d t + u^dk again. We have the equation for X 2 v 

L t X 2 i + [X 2 , □„.]« + X 2 N{\) + X 2 F = 0. (96) 

By Lemma 11, we can estimate the nonlinearity 

\\X 2 M(\)h* <e 2 + \\\Xi\ 2 + \i\ 2 + (|€| + lUr 1 )^! + \X 2 i\) + \X 2 a \(\i\ 2 + 1^)^ 
<e 2 + \\Xv\\ L ~\\Xv\\ L i + \\vf L4 + e\\X 2 v\\ Le + \\X 2 a \\m\\v\\m 
<e 2 + e\\Xv\\ H2 + \\v\\ 2 H1 +e\\X 2 n Hl 
< emax{e, ei} + e\\d 2 Xv\\ L 2. 

For the commutator, using (88), we can show that 

\\[X 2 ,L t ]i\\ L , < \\{X 2 ,a g ^\\ L2 + ||[X 2 ,yl( 7 )]£|| L2 + ||[X 2 ,^6]^|| L2 + \\[X 2 ,f^]v\\ L i 

< \\d(g*)d 2 Xv\\ L2 + \\d 2 (g*)d 2 n L2 + \\&(g e )dv\\v + \\d 2 (g*)v\\ L 2 
+ \\d 2 (f)di\\ L2 + \\X 2 f^v\\ L , + \\Xf^Xv\\ L2 

< e\\8 2 Xv\\ L2 + ||9 2 (.g e - ^)|| L6 |||9 2 5| + \dv\ + \v\\\ L s + \\d 3 (g* - h e )\\v\\dv\\ L oo + e 2 

< e\\d 2 Xi\\ L2 + \\d 2 r\\m (e + \\d 2 v\\ m ) + + ^ 
<e\\d 2 Xi\\ L2 +e 2 \\d 3 i\\ L2 +e 2 . 



Using the improved estimate (95), similarly to XF, we can estimate X 2 F as follows 



II^IU' < e 2 + \dh\ + \\d 2 g e d 2 <p s \\Li + W&g'dfaWv + \\d 2 g e (\dv\ + \v\)\\ L , + W&tfvh 



< e max{ e , ei } + ||c?V - h')\\ L , + \\d 2 (g* - h*)\\ H , 

< emax{e, ei}. 

Thus apply Lemma 12 to equation (96). We get 

\\dX 2 i\\ 2 L2 < max{e 2 , e 2 } + e' 1 f e 2 \\d 2 Xi\\ 2 L2 + e 4 ||9 3 5||| 2 ds < max{e 2 , e 2 } + e f \\d 3 i\\ 2 L2 ds. 

Jo Jo 

To retrieve the full estimates ||i9 3 i)||^ 2 , we apply Lemma 13 to the equation 

L e dv + [d, L € ]v + dJ\f(\) + dF = 0. 

We can show that 

\\dv\\ H 2 < max{e, ei} + \\d 2 XH\\ L2 + || [d, L e ]v + dAf(\) + dF\\ L2 

< max{e, Cl } + \\d 2 Xi\\ L2 + \\(\v\ + \v\p~ ^Qdf^ \ + \dv\)\\ L , 
<max{e, ei } + ||c> 2 X{S|| L2 . 

In particular, we have 

\\Xv\\h, <max{e, ei } + ||aX 2 ii|| L2 . 
Therefore from the estimates for dX 2 v we have obtained above, we can show that 

\\dX 2 i\\ 2 L2 < max{e 2 , e 2 } + e /* \\d 3 i\\ 2 L2 ds < max{e 2 , e 2 } + e f \\d 2 Xi\\ 2 L2 ds 

Jo Jo 

<max{e 2 ,e 2 } + e / \\dX 2 i || 2 2 ds, Vt < T/e. 
Jo 
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Then Gronwall's inequality implies that 

\\d 3 v\\ L 2 <max{e,ei} + ||aA- 2 t;|| ia <max{e,ei}, Vie [0,T/e\. 

This together with estimates (94) proves Proposition 7. 
Finally, the estimates (95) imply that 

||d s (0 - <t> s (x; \{t)))h* < \\d s (M*; Mm + e i@ W - cj> s (x; \{t)))\\ L , 
< max{e, ei} + |<9 2 7| < max{e, ei} 

for all \s\ < 3. Hence we have finished proving Proposition 3. 

7 Proof of the Main Theorem 

We use boostrap argument to prove the main Theorem 4. Using the Fermi coordinate system, we have 
shown the existence of solution <f> of equation (36) as well as its properties in Theorem 5 and Proposition 3 
under the assumption (37), which could be viewed as a boostrap assumption for the full reduced Einstein 
equations (33). We consider the equations of tp e = g e — h e to improve this boostrap assumption and thus 
to conclude the main Theorem 4. 

7.1 Estimates of the Metric g e 

Let (<7 e , (j)) be a solution of the system (33) with initial data satisfying conditions (30), (35) on the space 
([0, T/e] x M 3 , h e ). We have shown in Lemma 2 that the difference ip e = g e — h e satisfies the following 
hyperbolic system 

-{gT^ap^u + «V + <*V + SQ^ = 2J 2 (T„„ - UrT ■ «£„), 

where <5Q M „, SZ^, SP^ are given in (34). We show in this subsection that 

Proposition 9. If ip e = g £ — h e satisfies condition (37), then 

\\d s+1 r\\ L 4t)<5 2 , Vt<T/e,\s\<2. 

The key observation that allows S < e q , q > 1 is based on the fact that the energy momentum tensor 
T^virf 1 } splits into soliton part, which travels along the timelike geodesic (t, uot), and the error term which is 
small by Proposition 3. When doing energy estimate, we multiply the equations by Xip e = (<9 t +«oV x)ip e ■ 
By using integration by parts, we can pass the derivative X to the soliton part of T M „[0]. This, according 
to (88), allows us to prove Proposition 9 for all 5 < e q , q > 1. 

Proof. Since the initial data {<f>Q,<f>\) satisfy condition (30), we conclude according to Theorem 5 and 
Proposition 3 that <j> decomposes as (86) associated to the modified curve X(t) such that the remainder 
v satisfies the estimates 

\\d s v\\ L i(t)<e, V|s|<3, t<T/e. 
Using the modified decomposition (86), we can write 

r„„ - \trT ■ <£„ =< dph 8 v <t> > +V{<j>)g% = T* v + 
with the soliton part given by 

TjL =< d^ s (x;Mt)),d^ s (x;\(t)) > +V(<j> s & *(*)))<£,, 

= dpfMu, + a M e(A)^e(A)/ 2 o + v{f uo ) g i v . 
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The error term T^ v is small by Proposition 7. In fact, we can show that 

E \\9 s T^h<t)<e+ J2 ll^ +1 «IU 2 + ||a s (V(0)-V(0s))||L 2 + ||aV(V(0)-V(^))||L 2 

|s|<2 |a|<2 

< e + \\d 2 {9 e - h<)\\v + ||V"(0) - V"(0 S )|| L2 + ||V'(0) - V'(0 S )IU» 
<e+U-4> s \\L- <e, Vt<T/e. 

Here we recall that V(<j>) is given in line (5) and p > 2. 

Since X = dt + u^dk is timclike, apply formula (63) with Y = X, j3 = 1 to the above hyperbolic 
system for ty* commuting with the vector field d s . We obtain the energy estimates 

\\dd s r\\h(t) < \\dd s r\\h(0) + fiWd'g^d'rW^ + Wd^SP^+SZ^ + SQ^W^WXd'rW^ds 

Jo 

+ ef \\dd^r\\hds + S 2 f \\d s T*\\ L 4Xd s r\\L>ds + 5 2 f f d s T^ ■ XVp^dvol 
Jo Jo Jo Jr 3 

where |si| < \s\. Recall that h e (t,x) = h(et,ex). By assumptions (21), we conclude that 

HMd s+2 /i £ || L oo + \\\x\{dh e )(dh £ )\\ L ~ + ||a s+1 /i e || L =o < e, \^\ + |9V e | < e. 
By the definitions of SQ^, SZ^ U , 5P^ given in (34), for \s\ < 2, we can estimate 

\\d s sz^\\ L 2<e J2 \\d si+l rhi + J2 \\\x\d" +2 h e -\ x \- i r\\L'<t ]T w i+1 r\\v, 

Si<S — 1 S2<S Si<S— 1 

\\d s 6P^\\ L2 < e ]T H^+VIl' + \\\x\d s (dh* ■ dtf) ■ M-VIIl^ < e ]T ll«* 1+ V e |U», 

S±<S S\<S 

\\d s 6Q^\\ L 2 < e £ ||9 S1+ V1I^ + IIM0W • • M" 1 ^!^ + II^V ■ ^Vll^ 

Si <s 

<^E \\d si+1 r\\L* + \\d 2 r 

Si <S 



Here we use Lemma 5 to bound |||x| 1, c ||l 2 - 

For the soliton part, first notice that with the modified curve A, we can compute 

Xf UQ (z{x; A)) = -V z f U0 A UQ u, \u(t)\ < f \j\ds < e. 

Jo 

Therefore by Proposition 8 and inequality (95), we have 

\\Xd-\f^z)(l + \x\)\\ L2 (t)<\u\ + \j\ + \l\ + \dh\<e, Vt<T/e, \s\ < 4. 
Using integration by parts and Lemma 5, we can show that 

f f d s T^ ■ Xd s ^dvol 
Jo Jr 3 

< \\d s T*„ ■ d s % v \W{t) + \\dX» ■ 9 s %Mm + f UdXu ■ d s % v \Wds 

Jo 

< wa'T^i + \x\)\\»\\d>p^{i + m^im*) + \\d s T^(i + \x\)\\ L2 \\d°^Ai + Mr'um 

+ f \\Xd s T^(l + |x|)|| L2 ||d^(l + Ixir^ds 
Jo 

< \\d s+i r\\L>(t) + \\d s+i r\\L>(o) + e f \\d s+i rwd S . 

Jo 
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For \\d 2 g e • d 2 tp e \\ L 2, note that 

ii^y • d 2 r\\ L > < e 2 \\d 2 r\\L* + \\d 2 r ■ < t 2 \\d 2 r\\ L i + \\d 2 r\\ 2 m- 

Hence by the conditions (35), we can show that 

\s\<2 | S |<2 J ° J ° J ° 

+5 2 \\d s+i rh<t). 

Since we have assumed 

ll^+Vl^ <2e, Vt<T/e, 
Gronwall's inequality then implies that 

\\dd s r\\L*(t)<5 2 , Vi<T/e, |s|<2. 

Thus the proposition follows. □ 



7.2 Proof of Theorem 4 

The foliation of the spacetime is not preserved under the change of coordinate system constructed in 
Lemma 1. Since the argument of Theorem 5 relies on the new Fermi coordinate system, we first extend 
the given vacuum spacetime ([0, T] x R 3 , h) to ([0, T + 2Si] x M. 3 ,h) for some small positive constant Si. 
This can be obtained due to the assumptions (21) together with the local existence result for Einstein 
equations [4]. The metric h still satisfies condition (21) but with some new constant Kq. Therefore Lemma 
1 implies that one can choose a Fermi coordinate system (s,y) <G [0, cq{T + Si)] x K 3 on a subspace M 
such that 

M = [0,r]xK 3 cMc[0,T + 2i ]xI 3 , (97) 

where c = C(uh(0), h), depending only on the initial data Ao € A sta b(0)(see the definition in the proof of 
Lemma 1). We now can identify M with the space [0, Cq(T + Si)] x R 3 with the Fermi coordinate system 
(s, y). Then on the rescaled space [0,cq(T + Si)/e] x M 3 , the hyperbolic system (33) with initial data 
described in Lemma 2 admits a unique solution (g e ,(f>) on [0, t*) x R 3 for some small time t*. Moreover, 
Proposition 9 implies that under the assumption 

sup V \\d a r\\ L ^)<2e 2 , 

°^**|a|<2 

we in fact can show that 

sup V \\d a r\\ L 2(s) < C 3 S 2 < C 3 e 2 «(or C 3 e 2 e 2 if S = e e) 

for some constant C 3 independent of e. Note that q > 1. Additional to the requirement on e in Theorem 
5, if we choose e(or e ) such that 

C 3 e 2 «- 2 (or C 3 e 2 ) < 1, 

then we can improve the boostrap assumption (37). This also implies that the solution {g e ,(f>) of (33) 
can be extended to the whole space [0, cq(T + Si)/e\ x R 3 such that 

£ ||dV-^)|| L2 ( S )<<5 2 , Vs<c (T + Si)/e, 

1<I/3|<3 

ll^-0s(y;A e ( S )))|| i2 ( S )<e, Vs<c (T + Si)/e, 

|/3|<3 
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where we define the modulation curve X e (s) as follows 

A £ (s) = (w(es), e~ 1 d(es),e- 1 ((es) + w s, u a + u(esj) e A stab . 
Moreover, this curve is close to the given time like geodesic (s,uqs) in the sense that 

|C(s)| + |w(«) - w(0)| + < e, Vs < c (T + fc). 

We now use these results to construct solutions of the Einstein equations (1). Under the Fermi 
coordinate system, rescale the spacetime ([0,co(T + #i)/e] x R 3 ],g e ,<j)) to ([0,co(T + x R 3 ,g, C ) in 
the following way 

9(s,y) = g e (s/e,y/e), (j) e {s,y) = 5<t>(s/e,y/e). 

Making use of Lemma 2, we conclude that ([0, co(T + Si)] x R 3 ,g, <p e ) = (M,g,(j> e ) solves the Einstein 
equations (1) and satisfies the initial data (S , g, K, cj) , (f>\) given in Theorem 4. The Fermi coordinate 
system (s,y) on M leads to a foliation S T of M in a natural way S T := {s = t} such that for such 
foliation we have the following estimates for the solution (g 7 (f> e ) 

E ^W^(9-h)\\ L ^s)<S 2 , \/s<c (T + Si), 

1<|/3|<3 

E e a - f l|9 a (5-y-M2//e;A e ( S /e)))|U 2 ( S )<e, V S < c (T + h). 

|a|<3 

Let the C 1 curve 

A(s) = (u(s),6(s), C(s) + u s, u + u(s)) € A sta b 
be defined from A £ (s), which has been given above. Recall Definition 1 for 4> e s (y; A(s)). We note that 

s (2/A;A e ( S /e)) = ^(2/;A( S )) 

if h(s, C(s) + mqs) = too- However, since 

lC(s)|<e, h(s,u s) = mo, 

we conclude from (42) that 

E £Q_i II 9 " (^(2/5 AW) - M2//e; A e (s/e))) (s) < \h(s, C(s) + u s) - h(a, u s)\ < e 2 . 

|a|<3 

In particular, we have shown that 

||*-V(*,i/)-0s(y;A(*))||ffs +e\\S- 1 <f>l( S ,y)-rs(y^(s))\\H ? <e, V S <c (T + (5i). 

This proves estimates (25), (26). 

Finally, since the space M(diffcomorphic to [0,co(T + #i)]) can be viewed as an extension of Ai by 
(97), restricting the solution (M,g, <jf) to Ai, we obtain a solution (A4,g, 4> e ) of (1) as well as a C 1 curve 
A(s) = (uj(s),6(s), ((s) + uqs, uq+u(s)) such that (27) holds. This solution is unique up to diffcomorphism 
by a result in [9]. This completes the proof of the main Theorem 4. 
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8 Existence of Initial Data 



In this section, we discuss the existence of the initial data (R 3 , g, K, </> , <j>\ ) satisfying the conditions in 
Theorem 4. 

Let (U. 3 ,h,k) be the given initial data for the vacuum spacetime (A4,h), satisfying the vacuum 
constraint equations 

R{h) - \k\ 2 + (trk) 2 = 0, - V.trfc = 0, (98) 

where R denotes the scalar curvature on (R 3 ,/i), V is the covariant derivative with respect to h. Let 
{x\(xi,Xz,X3)} be a coordinate system on R 3 . Assume (j>Q, 4>\ are given functions on R 3 of the form 

4>l(x) = SMx/e) = 8M-/e), <t>\{x) = Se^M^M = <^~Vi(-/e), 

where 5 = e q , q > 1 or 5 = e$e. We want to show that there exists a Riemannian metric g and a symmetric 
two tensor K on R 3 satisfying the Einstein constraint equations 

R(-g) - \K\ 2 + (trK) 2 = b\-\\4>tf + |V^o| 2 + 2V(tfo))(-A). (QQ) 
V^Kij ~ VitrK = S 2 e- 2 < Vjfo > (-/c), 

as well as the estimates 

\\V{9-h)\\Hi(R3) + \\K-k\\ H i(R3) ^CifatuhtyPe- 1 . (100) 

Here V is the covariant derivative for the unknown metric g and the function V is defined in (5). 

We define the weighted Sobolev space H s,w on R 3 

H\\h..«:= ( £ / \d m <P\ 2 (l + \x\ 2 r +m dx) , 

\v< m <s J ^ 3 J 

and weighted Holder space C°' w 

|H| C o,„ :=sup{(l + Hn^)|}. 

X 

We define the metric space M s ' w on R 3 as follows: 

M s ' w :— {Riemannian metric g;gij — (toq)^- G H s ' w }, 

where too is the Euclidean metric on R 3 , that is, (mo)u = 1, (rao)ij = if i ^ j. We have the following 
existence result of the initial data (R 3 ,g, K 7 (f>Q , (f>\) . 

Theorem 6. Let (Eo,h, k) be the initial data for the vacuum spacetime (A4,h) such that h G M 4 ' -1 , 
k G H 3 ' a , trk — 0. Assume the matter field 0o € H 3 ' -1 , <f>\ G H 2 ' . Then there exists e > such that 
for all e < e ; there exists a Riemannian metric g and a symmetric two tensor K satisfying the constraint 
equations (99) and the estimates (100). 

Remark 9. The method here also applies to the case trk = constant which has been studies in [8]. 
However, the assumption k G H 3 ' together with trk — constant imply that trk — 0. For the general case 
trk 7^ 0, see the work of Corvino and Schoen [11]. 

The existence of g, K has been shown in [5], [8] by using implicit function theorem. For consistence, 
we repeat the prove. However, the difficulty here is to show that g, K obey the estimates (100) for all 
5 = e q 1 q>lor5 = e e, in particular the estimate 

\\V(g-h)\\ L2 + \\K-k\\ L2 <C5 2 e?. 
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The approach of previous works can only imply 

\W9-h)\\ L > + \\K-k\\ L 2 <C5 2 e-^. 

We improve this estimate by relying on a Hardy's inequality(see Lemma 16) for the first order linear 
operator JV (defined below), which does not have any nontrivial kernel in the class H 2 '^ 1 . Before proving 
this theorem, we make a convention that A < B means A < CB for some constant C depending on h, 
00, 0i- 

Remark 10. Suggested by our argument for the main theorem, the previous results in [5], [8] may imply 
the above existence theorem if we consider the rescaled constraint equations. However, the estimates 
depends on the iJ 4,_1 norm of the given metric h. After scaling, the i? 4 '" 1 norm of the scaled metric 
h e {x) = h(ex) depends on e. In fact, it can be shown that the scaled norm has size e~2+ w = e~2. Hence 
considering the scaled constraint equations can not lead to the above theorem directly. 

We denote the Hamiltonian constraint 

H((g,K), (0 O , 0i)) = R{9) - \K\ 2 + (trK) 2 - S 2 e' 2 d^ 2 + |V0 O | 2 - 2V(0 o ))(-/e), 
and the momentum constraint 

M((g, K), (0o, 0i)) = Vf. - VtrK - 5 2 e~ 2 < cf> u V0 O > (■/*)> 
where V is the covariant derivative with respect to g. Define the spaces 

X :={(g-m ,K)\geM 3 >-\ K E H 2 > }, 
Y := {(0o,0i)|0o eH 3 '~\ 4>i eff 2 - }, 
Z:={{p,J)\ P EH^\ JeH 1 ' 1 }, 

where p is scalar function, J is vector valued function on R 3 . We define the constraint map 

$ : X x Y -> Z, 

(g - m , K) x (0o, 0i ) (H((g, K), (0 O , <f>i)),M((g, K), (0 , 0i))). 

The fact that $ is a map from X x Y to Z follows from the multiplication and embedding properties of 
the weighted Sobolev spaces. We state Lemma 2.4 and Lemma 2.5 in [6] here. 

Lemma 14. We have 

3 3 

< C\\f\\ H s 1 ,™ 1 \\g\\ H s 2 .™2, s 1 +s 2 > s+-, w 1 +w 2 > w - -, 

||/Hco.«'<C||/||ff..», 8>\, w'<W+ i -, 

where the constant C depends only on s, w, s\, w\, s 2 , w 2 , w' . 
Let 

x = (h,k), y = (0,0), x = {g,K), y = <5 2 e - 2 (0 o , <£i)(7<0- 
Then the vacuum constraint equations (98) become $>(xo,yo) = 0. We define the linear map 

D<S>(x ,y ) :X^Z, 

(g,K)^(DH,DM) 

as the linearization of $>(x, y) at the point (x ,yo), which can be computed as follows 
DH = -A- h {tr- h g) + div- h (div- h g) - g ■ Ric(h) - 2k ■ K + 2(tr- h k)(tr n K), 
DM = VjKl - Vi(tr- h K) + ^ktytrzg - \w l V i9jl . 
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Here the covariant derivative V is for the metric h. The above formulae could be found in [5] or can be 
obtained by straightforward computations, noticing that the linearization of the connection is given by 

5T ij = 2 ( V 4 + V jff' - Vfltf) . 9i = 9ilh lj 

for symmetric two tensor and gij. Moreover, using Lemma 14, we can show that 

mijg + h, K + k), y ) - $((h, k), y ) - (DM, DM)\\z < II (g, K)f x . 

To apply the implicit function theorem, we must show that the linear map D$(xo,yo) is surjective from 
X to Z, that is, for any z = (p, J) G Z, the equations 

DU = p, DM = J (101) 

have at least one solution (g,K) G X. Notice that the above equations are underdetermined. The linear 
map D<fr(xo, t/o) has nontrivial kernel. We are instead looking for a solution of the form 

g=\\h, K = L v h-div{V)h-\\k (102) 

for some real function A and vector fields V on R 3 . Here Lyh is the deformation tensor of the vector 
fields V on (R 3 , h) defined as follows 

(L v h)ij = ViVj + VjVi. 

Note that trk = on R 3 . Using (98), the equations (101) are reduced to the following elliptic systems 

-A h \+\k\ 2 \ = 3k ■ L v h+^p, (103) 
div(L v h) = J. (104) 

We must show that the above elliptic systems have a unique solution (A, V) £ H s,w for any z = (p, J) E Z. 
Since the systems are splitting. We first consider the second equation (104) which is independent of A. 
The following lemma indicates that the operator div(L^h) is injective from H 2 '^ 1 to H ' 1 . 

Lemma 15. Let h G M 4 ' -1 . If div(L v h) = 0, V G H 2 ' 1 , then V = 0. 

This result has been proven in [7] . From a geometric point of view, a killing vector field V is uniquely 
determined by V\ p , VF| p . The condition div(Lvh) — 0, V G H 2 ^ 1 implies that V is killing and V, \7V 
vanish at infinity. Hence V vanishes everywhere. However, we give another proof inspired by the method 
in [10], see Theorem 3.3 there. 

Proof. Since V G if 2 ' -1 , we have 

= / h{div{L v h),V)da = -\ [ \L v h\ 2 da 

by approximating V with vectorfields V n G Cfi°. Hence Lyh = 0, that is, V is killing. In local coordinates, 
we have V^Vj + Va^-V^ = 0. In particular, we can compute 

V di V dj V k = -R(V,d i ,d j ,d k ), (105) 

where di is the vector field d Xi , R is the Riemann curvature tensor defined as follows 

R(di,di,dj,dk) = h(\7 dl \7 d .dj - V di V^dj, d k ). 

Using Lemma 14 by taking w' = j < — 1 + 2 + |, we have 

iRidiMM < \W 2 h\ < (1 + \x\)~i, \V\ < (1 + \x\)~i. 
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Hence by Lemma 5 (or Poincare inequality), we can show that 

||(1 + \x\)V 2 V\\ lhm3/Bro) < (l + R)-i\\(l + \x\)- l V\\ LHWL3/BRa) 

<{l + R)-H{l + \x\)V 2 V\\ L ^ /BRa) 

for any ball B Ba with radius R . Choose R large enough. We can conclude that V is vanishing outside 
the ball Br . 

Now consider the set 

S:={i£l 3 , V\ x =0, W|*=0}. 
We show that the set S is open. In fact, let x £ S. Notice that V £ C°. We have 

/ \V\ 2 {x + SLo)ds<r 2 f |W| 2 ds<r 4 / \W 2 V\(x + suj)d S < r 4 f \V\ 2 (x + soj)ds, Vcj G § 2 , 
Jo Jo Jo Jo 

where we have used the equation (105). Choosing r small enough, we can show that the ball B r (x) C S. 
Hence S is open. Notice that S is closed and nonempty. We conclude that V = on M 3 . □ 

This lemma also implies that the first order linear operator L^h is injective. We prove a Hardy's 
inequality for this operator, which will be used to improve estimates for ||(VA, VV)||l2. 

Lemma 16. Assume h G M 4 " 1 , V G H 2 ' 1 . Then 

IIN"V|| L 2 < \\L v h\\ L *. 

Proof. Choose Rq such that 

\Ric\ < ^(l + \x\)- 2 , \x\>R . 

We claim that 

\\V\\ LHBRa) <\\L v h\\ L ,. (106) 
In fact, if the above inequality does not hold, assume V n G H 2 ^ 1 such that 

l = \\V n \\ L 2 {BRo) >n\\L Vn h\\Li. 

Integration by parts, we have 

l\\L v h\\ 2 L2 = f VW^iVj + VW^iVjda 
* Jr 3 

V^ViVj - V^ViVj + V^i^Vj - V^i^Vjda 

= [ \VV\ 2 + \div(V)\ 2 - Ric(V,V)da. 
Jr 3 

Hence we conclude from the assumptions on V n that 

HVl^Hta < 1 + \\Ric\\ c o\\V n \\ L 2 {BRo) + -^||(1 + \x\)- l V n \\ L 2 {m? , /BRa)l 
which, by using Lemma 5, implies that 

||W„||i2 <2 + 2||ffic|| c o. 
Therefore, there exists vectorfields V such that (up to a subsequence) 

W„ W weakly, V n V strongly in L 2 (B ). 
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In particular, we have 

\\L v h\\ L 2 < Jim \\L Vn h\\ L 2 = 0, \\V\\ L 2 (Bao) = 1. 

That is V is killing vector field. Moreover, by (105), we have V G H 2 '^ 1 . Then Lemma 15 implies that 

V = 0, which contradicts to ||V"||i,2( Bji j = 1. Therefore the desired inequality (106) holds. Thus , for 

V G H 2 '^ 1 , we can show that 

\\\x\- 1 v\\ 2 L 2<\\vv\\ 2 L 2 = 2\\L v h\\ 2 L2 + f m c (v,v) - \div(v)\ 2 da 

< 2\\L V h\\ 2 L2 + iiffi C || C o||T/|| i2(BRo) + i||(i + wr'vw 

which implies that 

iiM-^iii* < mvhwh + n^c|| C o||y|u 2(BKo) < \\L V h\\ 2 L2 . 

□ 

Remark 11. An alternative approach for Lemma 15 and Lemma 16 is to use continuity argument. We 
sketch the prove here. Let 

L t (V) = L v (th + (l-t)m ), ie[0,l]. 

If for some t G [0, 1] such that 

WW-'VWv <C \\L to (V)\\ L 2, VVeH 1 '- 1 , 

then we can show that 

\\L U) {V)\\ L 2 <C(h)\\L t (V)\\ L 2 

for t close to to- Since Lemma 16 follows from Lemma 5 if h = tuq, we thus conclude that Lemma 16 
holds for all h € M ' , V G H 2 '^ 1 . In particular, Lemma 15 follows from Lemma 16. 

We now proceed to prove Theorem 6. By Lemma 15, the operator L^(th + (1 — t)mo) is injective 
from to II s ' 1 , for s — or 1. Now, for t = 0, the operator L^m — 2A is a diffeomorphism 

from iJ s + 2 > _1 to H*' 1 , see Theorem 5.1 in [6]. Hence the method of continuity [14] implies that the 
operator L^h is a diffeomorphism from H s + 2 >- 1 to H"' 1 . In particular, there exists a unique solution 

V G H s + 2 >~ { of (104) such that 

\ H s + 2,-l <\\J\\„.,1. 



Moreover, by Lemma 16, we can show that 

\\L v h\\ 2 L 2 = -2 f h(V,div(L v h))da < \\\x\- x V\\ L 2 ■ \\\x\J\\ L 2 < \\L v h\\ L 2\\\x\J\\ L 2, 
which implies that 

\\L v h\\ L 2<\\\x\J\\ L 2. (107) 

Having obtained V , for equation (103), by Theorem 6.6 in [6], the operator — V^ + |fc| 2 is a diffeomorphism 
from H s + 2 ^ 1 to II s ' 1 , for s — or 1. By Lemma 14, k ■ Lyh G H 1 ' 1 . Therefore, there exists a unique 
solution A of (103) such that 

||A|| ff . +2 .-i < \\k ■ L v h\\ H .,i + \\p\\ H .,i < \\p\\ H .,i + \\J\\ H ^, a = 0,l. 
Moreover, multiply equation (103) by A and integration by parts, we can show that 

||VA||| 2 + ||fcA||| 2 = - A (3k ■ L v h+ ^ da < M^X]^ (\\\x\ P \\ L 2 + \\L v h\\ L 2) , 



52 



where ||x|fc| < 1 by the assumption that k £ H 3 ' . Using (107), we get a similar estimate 

l|VA|| i2 < HIsHI^ + lllxlJll^. (108) 

Now we have shown that the linear map D$(xo,yo) is surjective from X to Z. Hence the Banach space 
X can be decomposed as X = X\ + X 2 such that D$>(xo, yo){X2) = 0, £)$(xo,j/o) is a diffeomorphism 
from X\ to Z. In particular, for y = 8 2 e~ 2 ((f>o, 0i)(-/e), the implicit function theorem shows that there 
is a solution 

(g,K) = (h + g,k + K)£X 

of (99) if e(or eo) is sufficiently small, depending only on h, cj>o, <f>i. Moreover, we can require that (g, K) 
is of the form (102) for (A, V) £ Z. Therefore 

11(^*011^.-1 < <5 2 e- 2 ||(|^i| 2 + |Wo| 2 - 2V(^ ))(-/e)|| jr o.. + || < cf> u Vfo > (-/e)|| ff o,i < 5 2 e^ , 
which implies that 

6 X -i (||V 2 (.g - h)\\„ + \\V(K - k)\\ Ll ) < e-*||(A, V)\\ H ^ < 5 2 e~\ 
We must remark here that in local coordinate {x} we have </>o(7 e ) = <Ao(x/e). Similarly 

II (A, V)\\ H 3,-i < <5 2 e - 2 ||(|^i| 2 + |Wo| 2 - 2V(^ ))(-A)|| ff M + II < <t>u Wo > (-A)llfltM < <S 2 ^, 
which shows that 

e 2 " 1 (l|V 3 (.g - + ||V 2 (if - fc)|| L2 ) < ei\\(\,V)\\ H3 ,-i < 6 2 e~\ 

It remains to estimate ||V(g — h)\\ L 2 + \\K — k\\ L 2, which we rely on (107), (108). Note that 

\Vh\ + \k\ < \x\-\ 

Making use of the multiplication properties of H 3 ' w in Lemma 14, we can show that 
\\V(g-h)\\ L 2 + \\K-k\\ L2 < \\VX\\ L 2 + \\L v h\\ L2 

< 5 2 e- 2 (|||x|(|^| 2 + |Wo| 2 - 2V(^o))(-/e)|U a + |||x| < <j> u V^ > (-A)lM + \W{g,K)\\ H ^ 

< S 2 d + WgW 2 ^ + \\K\\ 2 H1 , < S 2 d + {5 2 e^f < S 2 e^, 

where Af(g, K) is the nonlinear term in equation (99) forg — h + g,K~k + K. Hence 

\\V(g h)\\ H 2 + \\K- k\\ H 2 = ]T e ' Qhi (l|V Q+1 (ff - h)\\v + \\V a (K - fe)|U=.) < 5 2 e-\ 

|a|<2 

That is the solution (g,K) satisfies the constraint equations (99) as well as the estimate (100). 
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